MICROCOPY  RESOLUTION  TEST  CHART 

1»'  ON»i  SuACAU  Of  STAHOMm  -  -  * 


COORDINATED  SCIENCE  LABORATORY 

College  of  Engineering 
Decision  and  Control  Laboratory 


AD-A190  499 


REPRESENTATION, 
ANALYSIS,  AND 
DESIGN  OF 
MULTIRATE 
DISCRETE-TIME 
CONTROL  SYSTEMS 


Kevin  Lee  Buescher 


DT  JC 

ELECTE 

MAR  0  3 1988 

°sJD 


I 


UNIVERSITY  OF  ILLINOIS  AT  URBANA-CHAMPAIGN 


Approved  for  Public  Release.  Distribution  Unlimited. 


88  3  03  07  S 


l«TOn«.miiHr«Mr.CTP.iacitwjM4 


REPORT  DOCUMENTATION  PACE 


I«.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


WMiMJEZ. VJSlrk 


2b.  OECLASS1FICATION  /  DOWNGRADING  SCHEDULE 


3  DISTRIBUTION /AVAILARIUTY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


A.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

UILU-ENG-88-2211  (DC-101) 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Coordinated  Science  Lab 
University  of  Illinois 


6c  ADORESS  (City,  State,  and  2/P  Code) 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

Office  of  Naval  Research 


7b.  ADDRESS  (Oty,  State,  and  ZIP  Code) 


1101  W.  Springfield  Ave. 
Urbana,  IL  61801 


800  N.  Quincy  St. 
Arlington,  VA  22217 


8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION  joint  Services 
Electronics  Program 


Be  ADORESS  (Oty,  State,  and  ZIP  Cod*) 
800  N.  Quincy  St. 
Arlington,  VA  22217 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Of  applicable)  N00014-84-C-0149 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO. 


WORK  UNIT 
ACCESSION  NO 


it.  TITLE  (Inclod*  Security  Classification) 


REPRESENTATION,  ANALYSIS,  AND  DESIGN  OF  MULTIRATE 
DISCRETE-TIME  CONTROL  SYSTEMS 


12.  PERSONAL  AUTHOR(S) 


BUESCHER,  Kevin  Lee 


13a.  TYPE  OF  REPORT 

Technical 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM _  TO 


1*.  DATE  OF  REPORT  (Tear.  Month.  Day) 

1988  Feuruarv 


COSATI  COOES 


GROUP  SUB-GROUP 


1 8.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  idenbfy  by  block  number) 

multirate  discrete-time  systems,  control  systems; 


sr 


discrete-time  representation 


ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


Multi  rate  discrete-time  systems,  systems  in  which  sampling  and  d i sere te- 1 i me 
calculations  are  performed  at  two  or  more  rates,  arise  in  a  variety  of  appl i ca t i ons . 
The  multirate  character  of  the  -system  may  be  intrinsic,  due  to  digital  subsystems 
operating  at  multiple  rates.  Alternately,  the  multirate  nature  of  the  system  may 
be  induced  by  the  addition  of  sensors,  actuators,  and  d i screte-t ime  control 
structures  at  different  rates  for  the  purposes  of  economy  or  performance,  i 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

EJunclassifieoajnumited  □  same  as  rpt.  □  otic  USERS 

21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 

22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

22b.  TELEPHONE  (Include  Area  Code) 

22c.  OFFICE  SYMBOL 

OO  FORM  1473, 84  MAR 


L 


83  APR  *<111100  may  0*  used  until  axhaustad. 
All  othar  aditiont  are  obsoltt*. 


SECURITY  CLASSIFICATION  OF  THIS  PACE 

UNCLASSIFIED 


jnunik’xvxwTi 


REPRESENTATION,  ANALYSIS,  AND  DESIGN  OF 
MULTIRATE  DISCRETE-TIME  CONTROL  SYSTEMS 


BY 

KEVIN  LEE  BUESCHER 
B.S.,  University  of  Illinois,  1986 


THESIS 

Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Electrical  Engineering 
in  the  Graduate  College  of  the 
University  of  Illinois  at  Urbana-Champaign,  1988 


Urbana,  Illinois 


ACKNOWLEDGMENTS 


The  author  of  this  thesis  is  most  appreciative  for  the  assistance  of  his  academic  mentor  and 
his  cohort  in  sowing  the  seeds  of  multirate  dissent.  Professor  J.  W.  Grizzle.  The  author  would 
also  like  to  thank  D.  Connors,  S.  Cusumano,  T.  Ting,  and  S.  Zabin  for  their  timely  advice  and 
helpful  comments.  Finally,  the  author  is  grateful  for  the  support,  monetary  and  otherwise,  of  the 
Decision  and  Control  Laboratory  of  the  Coordinated  Science  Laboratory,  Ford  Motor  Company, 
and  his  parents. 


m 


TABLE  OF  CONTENTS 


CHAPTER  PAGE 

1.  INTRODUCTION  . 1 

2 .  TIME-INVARIANT  REPRESENTATION  OF  MULTIRATE  SYSTEMS . 3 

2.1.  Systems  Admissible  for  Analysis  . 3 

2.2.  The  Least  Common  Period  of  a  Multirate  System . 6 

2.3.  Expansion  of  Discrete-Time  Signals  and  Systems  . 8 

2.4.  Expanding  Multirate  Sampled  Continuous-Time  Systems . 13 

2.5.  Connection  Matrices  . 18 

2.6.  The  T-expanded  Representation  of  a  Multirate  System . 21 

2.7.  Inclusion  of  Periodic  Subsystems  . 22 

2.8.  Modular  Expansion  of  Multirate  Systems  . 26 

3.  PERIODIC  REPRESENTATION  OF  MULTIRATE  SYSTEMS  . 29 

3.1.  The  t-varying  Representation  . 29 

3.2.  The  mx-varving  Representation  . 34 

3.3.  Corresponding  Representations  . 36 

4.  ANALYSIS  OF  MULTIRATE  SYSTEMS  . 42 

4.1.  Stability  . 42 

4.2.  Controllability  and  Reconstructibility  . 44 

4.3.  Duality,  Stabilizability,  and  Detectability  . 46 

5.  DESIGN  OF  MULTIRATE  CONTROLLERS  . 52 

5.1.  T-expanded  State  Feedback  and  Observers  . 52 

5.2.  M-varying  State  Feedback  and  Observers  . 55 

5.3.  Transfer  Function  Controller  Designs  . 58 

5.4.  Comments  . 59 

6 .  THE  VARIABLE  COMPONENT  METHOD  APPLIED  TO 

MULTIRATE  DISCRETE-TIME  SYSTEMS  . 60 

6.1.  A  Time-invariant  Result  . 61 

6.2.  Interpretation  of  the  Time-invariant  Result  . 63 


tv 


I 


6.3.  Application  of  the  Variable  Component  Method 


CONCLUSIONS 


APPENDIX  SELECTED  PROOFS 


REFERENCES 


SYMBOL  TABLE 


Symbol 

N 

Q 

Rnxm 

Z+ 

V 

n  mod  p 

SISO 

MIMO 

BIBO 

GCD 

LCM 

LCP 

T 

x 

P 

A’ 

Ay 

In 

Om 

Omxn 

♦  ♦  ♦ 


Meaning 

the  natural  numbers:  {1,  2,  3,...  } 
the  rational  numbers 
an  n  x  m  matrix  with  real-valued  entries 
the  nonnegative  integers:  {0,  1,  2,...  } 
for  all 

n  modulo  p:  p  x  (fractional  portion  of  n/p) 

Single-Input-Single-Output 

Multiple-Input-Multiple-Output 

Bounded-Input-Bounded-Output 

Greatest  Common  Divisor 

Least  Common  Multiple 

Least  Common  Period 

the  LCP  of  a  multirate  system 

the  fundamental  period  of  a  multirate  system 

T/x 

transpose  of  the  matrix  A 

i,j  ^  element  of  the  matrix  A 

the  n  x  n  identity  matrix 

an  m  x  m  matrix  of  zeroes 

an  m  x  n  matrix  of  zeroes 

end  of  theorem,  procedure,  example,  etc. 


CHAPTER  1 


INTRODUCTION 


Multirate  discrete-time  systems,  systems  in  which  sampling  and  discrete-time  calculations 
are  performed  at  two  or  more  rates,  arise  in  a  variety  of  applications.  The  multirate  character  of  the 
system  may  be  intrinsic,  due  to  digital  subsystems  operating  at  multiple  rates.  Alternately,  the 
multirate  nature  of  the  system  may  be  induced  by  the  addition  of  sensors,  actuators,  and  discrete¬ 
time  control  structures  at  different  rates  for  the  purposes  of  economy  or  performance.  Two 
practical  examples  of  designs  resulting  in  a  multirate  system  are  an  idle  speed  control  for  an  internal 
combustion  engine  (Powell  et  al.,  1987)  and  a  quadruplex  videotape  recorder  (Rao,  1979). 

Previous  studies  of  multirate  systems  began  with  the  study  of  multirate  sampled-data 
systems  by  Kranc  (1957),  which  employed  transfer  function  techniques  later  improved  by  Coffey 
and  Williams  (1966)  and  Boykin  and  Frazier  (1975).  These  works  presented  complicated  methods 
which  served  the  purpose  of  assessing  the  stability  of  a  system  composed  of  continuous-time 
subsystems  and  samplers  at  multiple  rates.  Kalman  and  Bertram  (1959)  used  state  space  methods 
to  study  sampling  systems  of  great  generality.  Although  Kalman  and  Bertram  demonstrated  that  a 
state  space  representation  exists  for  almost  any  hybrid  system  incorporating  a  variety  of  sampling 
schemes,  the  generality  of  the  systems  involved  precluded  the  presentation  of  systematic  and 
efficient  methods  of  obtaining  a  state  space  representation  for  a  given,  fixed  system.  Meyer  and 
Burrus  (1975)  introduced  the  concept  of  block  processing  to  perform  time  and  frequency  domain 
analyses  on  individual  single-rate,  periodic  time-varying  and  multirate  (inputs  at  one  rate,  outputs 
at  another)  digital  filters.  Recently,  Araki  and  Yamamoto  (1986)  applied  a  mutation  of  block 
processing  to  the  analysis  of  a  continuous-time  system  with  outputs  sampled  at  multiple  rates  and 
fed  back  through  a  constant  gain  matrix  to  the  inputs,  which  were  sampled  and  held  at  multiple 
rates.  Currently,  efficient  techniques  of  analysis  and  design  algorithms  for  nmltirate  systems 
composed  of  both  continuous-time  and  discrete-time  subsystems  are  not  available. 

Inspired  by  Kalman  and  Bertram  (1959)  and  beginning  with  methods  akin  to  the  block 
processing  used  by  Meyer  and  Burrus  (1975),  this  work  presents  a  systematic,  efficient  method  of 
obtaining  a  “time-invariant”  representation,  referred  to  as  the  T-expanded  representation,  of  a 
class  of  composite  multirate  systems  containing  both  continuous-time  and  discrete-time 
subsystems.  In  recent  years,  results  have  appeared  in  the  literature  pertaining  to  periodic  discrete- 


time  systems,  with  multirate  systems  cited  as  an  (unqualified)  example  of  such  a  system.  In 
response,  this  work  proposes  a  procedure  for  obtaining  a  periodic  discrete-time  representation, 
termed  the  M-varying  representation,  of  a  multirate  system.  Analysis  of  multirate  systems  via 
the  T-expanded  representation  by  time-invariant  techniques  is  shown  to  reveal  important 
characteristics  of  a  certain  M-varying  representation  of  the  multirate  system.  Examples  of  the 
design  of  controllers  for  multirate  systems  by  conventional  analytical  techniques,  with  minor 
modifications,  are  presented.  The  variable  component  method  is  applied  to  multirate  systems  to 
provide  an  iterative  means  of  building  composite  controllers. 

Chapter  2  develops  a  notation  suited  to  T-expanded  representations  and  details  procedures 
•seful  in  obtaining  a  T-expanded  representation.  The  notation  introduced  in  Chapter  2  is  used 
throughout  the  remaining  chapters  and  provides  a  means  of  discussing  a  variety  of  concepts  related 
to  multirate  systems  in  a  concise  manner.  In  addition  to  describing  the  behavior  of  a  multirate 
system  over  lengths  of  time  other  than  its  period,  the  periodic  representations  developed  in 
Chapter  3  provide  a  theoretical  tool  for  deducing  the  interperiod  behavior  of  a  multirate  system 
from  its  T-expanded  representation.  Chapter  4  discusses  the  stability,  controllability, 
reconstructibility,  stabilizability,  and  detectability  of  multirate  systems  and  emphasizes  the 
properties  which  a  periodic  representation  of  a  multirate  system  inherits  from  the  T-expanded 
representation  of  that  system.  Chapter  5  applies  the  developments  of  Chapters  2,  3,  and  4  to  the 
analytical  design  of  multirate  controllers.  A  time-invariant  technique,  the  variable  component 
method,  is  applied  to  multirate  systems  in  Chapter  6  via  the  T-expanded  representation.  The 
material  in  Chapter  5  and  Chapter  6  serves  as  an  example  of  the  manner  in  which  results  for 
standard  discrete-time  systems  can  be  extended  to  multirate  systems. 


CHAPTER  2 


TIME-INVARIANT  REPRESENTATION  OF  MULTIRATE  SYSTEMS 

This  chapter  presents  a  systematic  method  for  obtaining  a  single  time-invariant  state  space 
or  transfer  function  representation,  called  a  T-expanded  representation,  for  members  of  a  class  of 
composite  systems  which  employ  sampling  at  multiple  rates.  This  method  first  converts  the  block 
diagram  of  a  multirate  system  into  a  structurally  similar  block  diagram  of  time-invariant  discrete¬ 
time  subsystems  through  state  space  calculations  easily  performed  using  any  software  package 
capable  of  discretization  of  continuous-time  systems  at  a  single  rate  and  matrix  multiplication, 
addition,  and  composition.  The  reduction  of  this  single-rate  block  diagram  to  a  state  variable 
description  or  transfer  function  matrix  is  then  a  standard  problem. 

2.1  Systems  Admissible  for  Analysis 

The  methods  to  be  developed  will  apply  to  a  general  class  of  systems  with  multiple 
sampling  rates.  To  simplify  the  exposition,  the  multirate  systems  considered,  unless  otherwise 
noted,  will  be  as  follows: 

Assumption  2.1: 

a.  The  ratio  of  any  two  sampling  periods  in  the  system  is  a  rational  number. 

b.  All  samplers  are  in  synchronism  at  t  =  0  and  are  followed  by  zero-order 
holds  of  the  same  period 

c.  Discrete-time  subsystems  are  envisioned  as  sampling  their  inputs  and 
producing  zero-order  hold-type  outputs  (in  synchronism  with  the  samplers  in  part 
(b)  above  at  t  =  0.) 

d.  The  outputs  of  zero-order  holds  as  functions  of  time  are  continuous  from 
the  right  at  their  associated  sampling  instants. 

e.  The  entire  system  can  be  partitioned  into  subsystems  which  are  either 
discrete-time  or  continuous-time  subsystems  with  only  samplers,  summers,  and 
branch  points  (a  point  on  the  block  diagram  where  a  line  branches  into  two  or  more 
lines)  connecting  these  subsystems  together  on  the  block  diagram.  In  addition,  if 
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there  is  a  path  on  the  block  diagram  between  any  output  of  one  of  these  subsystems 
and  any  input  to  one  of  these  subsystems,  then  this  path  must  satisfy  one  or  more 
of  the  following:  the  path  originates  at  a  discrete-time  subsystem;  the  path 
terminates  at  a  discrete-time  subsystem;  the  path  passes  through  one  or  more 
samplers. 

f.  There  is  no  path  from  (to)  any  external  input  (output)  of  the  entire 
multirate  system  to  (from)  any  of  the  continuous-time  subsystems  in  pan  (e)  above 
which  does  not  pass  through  a  sampler.  In  addition,  there  is  no  direct  path  from  an 
external  input  to  an  external  output  of  the  multirate  system  which  does  not  pass 
through  a  sampler.  These  properties  will  be  loosely  stated  in  the  following  as  “the 
external  inputs  and  outputs  of  the  multirate  system  are  sampled.” 

g.  The  subsystems  in  part  (e)  above  are  linear,  time  invariant,  and  causal, 
and  may  be  MIMO. 

♦  ♦  ♦ 

Assumption  2.1(a)  is  the  necessary  and  sufficient  condition  for  a  multirate  system 
composed  of  time-invariant  subsystems  to  have  a  finite  period.  Satisfaction  of  Assumption  2.1(e) 
is  possible  for  essentially  any  connection  of  continuous-time  and  discrete-time  elements  and  may 
entail  mine:  block  diagram  manipulations,  to  be  detailed  later  in  this  section.  If  a  system  does  not 
satisfy  part  (f),  the  introduction  of  samplers  into  the  block  diagram  by  modelling  judgments  may 
serve  to  approximate  the  actual  system  adequately  and  satisfy  Assumption  2.1(0.  Assumption 
2.1(g)  will  eventually  be  relaxed  to  include  time-varying  periodic  discrete-time  subsystems.  In 
summary,  the  essential  assumptions  to  keep  in  mind  are  that  all  external  signals  (inputs/outputs)  are 
sampled  and  held  and  that  the  ratio  of  any  two  sampling  rates  in  the  system  is  a  rational  number. 

Although  the  term  “linear  multirate  discrete  composite  system”  more  aptly  describes  a 
system  satisfying  Assumption  2.1,  multirate  system  will  denote  such  a  system  hereafter. 

Example  2.1: 

The  multirate  system  shown  in  Figure  2.1  appears  to  violate  Assumption 
2.1(e)  due  to  the  direct  link  between  the  two  continuous-time  subsystems  O;  and 
O:;  yi(s)  passes  through  a  summer  and  a  branch  point,  but  not  through  a  sampler, 
before  reaching  U3(s).  Let  Oj  and  02  be  represented  by  the  transfer  function 
relationships 

yi(s)l  fNi(s)/D,(s)l  [U3(s)l 

=  U[(s)  and  y4(s)  =  [N'3(s)/D;(s)  N4(s)/D:(s)]  ;. 

y;(s)J  LNt2(s)/D[(s)J  Lu4(s 
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Figure  2.1.  A  multirate  system. 
Using  U3(s)  =  y i(s)  +  U2(s)  and  y3(s)  =  U3(s),  it  follows  that 


~y2(s)~ 

N2(s)/Di(s) 

0 

0 

-ui(s)- 

y3(s) 

= 

N^sVDKs) 

1 

0 

u:(s) 

_y4(s)_ 

_  N3(s)N1(s)/D2(s)Di(s) 

N3(s)/D2(s) 

N4(s)/D2(s)  _ 

_u4(s)_ 

Thus,  the  multirate  system  shown  in  Figure  2.1  may  be  partitioned  by  considering 
the  area  inside  the  dashed  line  in  Figure  2.1  as  a  single  continuous-time  subsystem, 
0 3,  as  shown  in  Figure  2.2.  The  multirate  system  in  Figure  2.2  satisfies 
Assumption  2.1(e). 


Figure  2.2.  The  system  in  Figure  2.1  redrawn. 


Example  2.1  provides  a  clue  to  a  technique  for  partitioning  a  multirate  system  so  that 
Assumption  2.1(e)  is  satisfied.  If  an  initial  partitioning  of  the  system  into  continuous-time  and 
discrete-time  subsystems  (connected  by  only  summers,  samplers,  and  branch  points)  does  not 
satisfy  Assumption  2.1(e),  then  by  elimination  a  path  must  be  connecting  a  continuous-time 
subsystem,  Cj,  to  a  continuous-time  subsystem,  02.  which  may  pass  through  summers  and  branch 
points,  but  does  not  pass  through  any  samplers  (0i  may  be  the  same  subsystem  as  03.)  Define  a 
new  continuous-time  subsystem,  03,  consisting  of  the  subsystems  0[  and  02  and  their  individual 


inputs  and  outputs  (u1.u4.y2.  and  y4  in  Example  2.1,)  except  for  the  offending  pair  in  question  (y  1 
and  U3  in  Example  2.1.)  The  subsystem  03  must  also  be  provided  with  additional  inputs  (U2  in 
Example  2.1)  representing  other  inputs  to  any  summers  the  path  passes  through  and  additional 
outputs  (y3  in  Example  2.1)  representing  other  outputs  of  any  branch  points  the  path  passes 
through.  This  process  can  be  repeated  until  the  entire  multirate  system  satisfies  Assumption 
2.1(e). 


2.2  The  Least  Common  Period  of  a  Multirate  System 


Given  a  set  of  N  nonzero  sampling  periods  {Ti,...,Tn},  assume  that  the  ratio  of  any  two 
of  these  sampling  periods  is  a  rational  number.  Consider  the  ratio  of  each  sampling  period  to  a 
particular  sampling  period,  say  Tj: 

— =  —  i  =  1  N 
Tj  qn  ’ 

for  some  ru,  qn  €  N,  where  the  Greatest  Common  Divisor  (see  Niven,  1980)  of  r^  and  qn  is  1 

for  each  i.  Let  Ri  be  the  Least  Common  Multiple  of  the  set  {r^ :  i  =  1,...,N},  X\  =  Tj/Ri,  and 

_  .  _  qiiRl  :  _  1  V 
Pli  _  >  t  —  1,...,N. 

1  ll 

Note  that  pn  e  N  because  Ri/ru  e  N.  Then 

The  process  of  obtaining  ti  and  the  pn's  will  be  referred  to  as  normalization  with  respect  to  Tj. 
Theorem  2.1: 

If  the  sampling  periods  are  normalized  with  respect  to  Tk  (giving  Tj  =  pkjtk, 
i  =  1,...,N)  and  again  with  respect  to  Tj  (giving  T;  =  pj;XJ?  i  =  1,...,N,)  then  Xk  =  Xj 
and  pki  =  Pji,  i  =  1,...,N.  (See  proof  in  the  Appendix.) 

♦  ♦  ♦ 

Thus,  the  sampling  periods  can  be  uniquely  represented  as  Tj  =  pjX,  i  =  1,—.N,  regardless  of  the 
particular  sampling  period  used  in  the  normalization  process.  Each  sampling  period  in  the  system 
is  a  multiple  of  x;  x  will  be  referred  to  as  the  fundamental  period  of  the  multirate  system. 

Since  the  subsystems  in  the  multirate  system  are  time  invariant,  the  periodicity  of  the  entire 
system  will  be  established  if  it  can  be  shown  that  the  samplers,  which  are  in  synchronism  at  t  =  0, 
are  again  in  synchronism  at  some  time  t  =  T.  For  a  system  with  sampling  periods  {Ti,...,Tn},  let 
P  =  LCM  { pi :  i  =  1,...,N).  At  times  t  =  mPx,  m  €  Z+,  all  samplers  in  the  system  will  sample: 

mPx  mP 


Tj 


■e  Z*  for  i  =  l.—.N. 


‘■'vlcv 


i 


4 


I 

5 

I 

\ 

i 

I 

l 

i 

1 

i 

i 

I 


Thus,  the  multirate  system  is  periodic  with  period  T  =  Pt;  in  fact,  this  is  the  Least  Common  Period 
(LCP)  of  the  multirate  system. 

Theorem  2.2: 

T  =  t  (LCM  {p,  :  i  =  1,...,N})  =  Pt  is  the  shortest  length  of  time  over 
which  a  multirate  system  with  sampling  periods  {T1,...,Tn)  is  periodic.  (See 
proof  in  the  Appendix.) 

♦  ♦  ♦ 

In  the  remainder  of  this  work,  the  symbols  T,  P,  T,  and  p*  will  be  implicitly  associated  with  the 
meanings  assigned  to  them  above.  A  procedure  is  now  given  for  normalizing  a  set  of  sampling 
periods  {Ti,...,Tn}  as  T,  =  p,x  and  finding  their  LCP. 

Procedure  2.1: 

a.  Arbitrarily  choose  a  sampling  period  Tk  from  the  set  and  express  the 

sampling  periods  as  T,  =  Ticq^i/rtc,  ,  i  =  1 . N  ,  where  q^r^,  e  N  and 

GCD  {qk„rki}  =  1  for  each  i. 

b.  Set  Rk  =  LCM  {ry  :  i  =  1 . N} 

c.  Set  x  =  T|c/Rk  and  pi  =  q^R^/ric,),  i  =  1 . N. 

d.  Set  P  =  LCM  { pi :  i  =  1,...,N)  and  T  =  Pt. 

♦  ♦  ♦ 

The  following  result  will  be  the  key,  in  Section  2.8,  for  establishing  a  certain  modularity 
property  of  multirate  systems. 

Theorem  2.3: 

Let  the  LCP  associated  with  {Ti,...,Tk)  be  T(k)  and  that  associated  with 
{ Ti,...,Tk,Tk+i }  be  T(k+1).  Then  T(k+1)/T(k)  €  N.  (See  proof  in  the 
Appendix.) 

♦  ♦  ♦ 

In  words,  the  effect  of  adding  a  new  sampling  period  to  a  multirate  system  is  to  increase  the  LCP 
by  an  integer  multiple. 

Example  2.2: 

Suppose  N  =  4  and  the  sampling  periods  are:  Ti  =  3/5  sec,  T;  =  1/7  sec, 

T3  =  1/3  sec,  and  T4  =  300  msec.  If  the  periods  are  normalized  with  respect  to  Tj, 
Procedure  2.1  yields 

Ti  =  (l/l)Ti,  T2  =  (5/2 1  )T  1 ,  T3  =  (5/9)Ti,  T4  =  (1/2)T; 


Pi  =  126(1/1)  =  126,  p2  =  126(5/21)  =  30,  p3  =  126(5/9)  =  70,  p4  =  126(1/2)  =  63 
P  =  LCM  { 126,30,70,63}  =  126  x5  x7  =  4410 
T  =  4410(1/210)  =  21  sec 

The  reader  may  wish  to  verify  that  pi,...,p4,t,P,  and  T  are  invariant  under 
normalization  with  respect  to  T2.T3,  or  T4. 

♦  ♦  ♦ 


2.3  Expansion  of  Discrete-Time  Signals  and  Systems 

The  members  of  the  class  of  multirate  systems  under  consideration  are  time  varying,  but 
periodic  with  period  T.  Knowledge  of  the  state  of  such  a  multirate  system  at  time  t  =  0  and  the 
inputs  to  the  system  over  the  interval  of  time  [0,T)  is  sufficient  to  determine  its  state  at  time  t  =  T 
and  its  outputs  over  the  interval  [0,T).  By  the  periodicity  of  such  a  multirate  system,  the  manner  in 
which  its  state  at  time  t  =  nT  and  its  output  over  [(n-l)T,nT)  are  determined  from  its  state  at  time 
t  =  (n-l)T  and  input  over  [(n-l)T,nT)  is  identical  for  each  n  e  N.  Thus,  the  behavior  of  such  a 
multirate  system  over  all  time  is  time  invariant  in  terms  of  the  description  of  the  system's  behavior 
over  one  period.  The  price  paid  for  this  time-invariant  description  is  that  all  input  values  to  the 
system  and  all  output  values  from  the  system  during  one  period  must  be  accounted  for.  This 
entails  expanding  a  single  input  line  into  many  fictitious  inputs  which  represent  its  values  over  the 
period. 

The  result  of  collecting  all  values  of  a  discrete-time  signal  over  amounts  of  time  of  length  L 
will  be  referred  to  as  the  L-expanded  version  of  that  signal.  The  signals  in  question  may  be 
vector  signals  with  individual  components  at  different  rates.  To  prevent  the  number  of  components 
in  the  L-expanded  version  of  a  signal  from  changing  with  time,  L  must  be  an  integer  multiple  of  the 
period  of  each  component  of  the  signal.  Discrete-time  signals  will  be  denoted  by  lower-case  letters 
and  the  same  letter  capitalized  will  represent  their  expanded  versions. 

At  this  point  it  is  convenient  to  introduce  the  concept  of  a  signal  bundle.  A  signal  bundle 
is  a  vector  of  discrete-time  signals  at  the  same  rate  which  may  be  a  portion  of  a  larger  vector  of 
signals  but  is  distinguished  in  some  manner  from  the  rest  of  the  vector.  The  block  diagram  of  the 
system  primarily  determines  the  grouping  of  signals  into  bundles.  Before  applying  the  methods 
presented  here,  certain  manipulations  will  be  performed  on  the  block  diagram  of  the  system.  The 
lines  drawn  as  inputs  or  outputs  of  subsystems  in  this  modified  block  diagram  will  each  be 
designated  as  a  bundle  of  signals;  each  line  may  actually  represent  an  entire  vector  of  signals.  Later 
developments  will  reveal  the  utility  of  using  bundles  of  signals  and  properly  explain  which  signals 
to  place  in  bundles. 
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Let  L/Tj  =  lj  e  N.  The  L-expanded  version  of  a  bundle  of  signals  u(nT]),  n  =  0,1,. 
the  vector  of  li  blocks  U(kL): 

"  uCkUTi)  - 

.  “((kh+DTi)  .  _ 


U(kL)  = 


,  k  =  0,1.... 


Lu^kh+li-DTjlJ 

U(kL)  will  be  thought  of  as  a  bundle  of  signals  with  rate  1/L.  The  set  (U(kL)  :  k  =  0,1,...} 
contains  all  of  the  values  in  {u(nT0  :  n  =  0,1,...},  so  the  expanded  version  retains  ali  the 
information  in  the  original  signal.  Capital  letters  will  denote  expanded  versions  of  signals 
hereafter,  and  the  capitalized  portion  of  an  expanded  signal's  argument  (L  in  this  case,  but  T  and  M 
later  on)  will  be  a  real  number  which  is  an  integer  multiple  of  x  and  denotes  the  interval  of  time 
over  which  expansion  has  been  performed.  As  a  slight  abuse  of  notation,  let  u(kL)  =  u(kliTi). 
Then  the  L-expanded  version  of  u(nTi)  can  be  written  more  conveniendy  as 

u(kL) 

u(kL+Ti) 

U(kL)  = 

_u(kL+L-Ti)_ 

The  expanded  version  of  a  signal  may  be  expanded  again  to  yield  an  expanded  signal.  The 
procedure  for  expanding  an  expanded  signal  parallels  that  for  expanding  a  normal  signal,  with 
expanded  bundles  treated  as  bundles.  As  an  example,  if  M/L  =  m  e  N,  then  the  M-expanded 
version  of  U(kL)  coincides  with  the  M-expanded  version  of  u(nTt)  and  is  given  by 

r  U(jmL)  -i  r  uijmliTO 

_ _  U((jm+l)L)  uCOmli+iyro  .  „ 


UfjM)  = 


,  j  =  0,1,. 


j_U((jm+m-l)L)J  [_u((jmli-t-mLli- l)Tj)  J 

To  see  the  equivalence,  note  that  the  first  component  of  U(jmL)  is  u((jm)liTi)  =  u(jml]Ti)  and  the 
last  component  of  U((jm+m-l)L)  is  u(((jm+m-l)li+li-l)Ti)  =  u((jmli+mli-l)Ti).  The 
equivalence  of  the  M-expanded  version  of  a  signal  and  the  M-expanded  version  of  the  L-expanded 
version  of  that  signal  obviates  the  need  for  additional  notation  to  distinguish  between  the  two. 

Let  y  be  a  vector  of  signals 

rxn 


where  each  Vj  is  a  bundle  of  signals  with  period  T;.  If  L  is  an  integer  multiple  of  each  T,  (  for 


example,  L  =  LCP  {Ti,...,Tq},)  then  the  L-expanded  version  of  y  is 

rYi(kL)- 

Y2(kL) 

Y(kL)  =  _  ,  k  =  0,l,...  , 

_Yq(kL)_ 

where  each  Yi(kL)  is  the  L-expanded  version  of  yj.  Although  Y(kL)  consists  of  many  of  the 
bundles  yi,  Y(kL)  will  be  thought  of  as  being  composed  of  the  q  bundles  Y;(kL).  Note  that  each 
bundle  of  signals  has  been  expanded  in  place.  Thus,  the  boundaries  between  bundles  have  been 
maintained. 

Repeated  iteration  of  the  equations  governing  a  simple  discrete-time  system  reveals  the 
behavior  of  such  a  system  over  units  of  time  of  length  L  and  forms  a  basis  for  determining  the 
behavior  of  more  complex  systems  over  units  of  time  of  length  L.  Consider  the  single-rate 
subsystem  with  realization  (A,B,C,D)  and  period  T;,  where  u(nT,)  and  y(nT;)  are  each  a  single 
bundle  of  signals: 

x((n+l)T,)  =  A  x(nT;)  +  B  u(nT{) 

,  n  =  0, 1 .  (2.1) 

y(nT0  =  C  x(nT;)  +  D  u(nT,) 

Since  iteration  of  a  discrete-time  state  equation  can  only  be  performed  an  integer  number  of  times, 
let  L  be  chosen  so  that  L/Ti  =  li  €  N.  Note  that  xfkL+ljT,)  =  x((k+l)L).  Expressions  for 
x((k+l)L)  and  Y(kL)  in  terms  of  x(kL)  and  U(kL)  are  desired.  By  direct  computation, 

x(kL+Tj)  =  A  x(kL)  +  B  u(kL) 
y(kL)  =  C  x(kL)  +  D  u(kL) 

x(kL+2T;)  =  A  x(kL+Tj)  +  B  u(kL+T,) 

=  A2  x(kL)  +  AB  u(kL)  +  B  u(kL+T;) 

vdcL+Ti)  =  C  x(kL+Tj)  +  D  u(kL+Ti) 

=  CA  x(kL)  +  CB  u(kL)  +  D  u(kL+T,) 

l. 

x((k+l)L)  =  x(kL+l,Ti)  =  A1'  x(kL)  +  X  [Aa''m)B  ufkL+lm-DT,)] 

m=l 

!.-l 

y(kL+(li-l)T1)  =  CAGi'1)x(kL)-t-Du(kL+(!i-l)Ti)+  X  CAn'  m  l)B  u(kL+(m-  1)T,)  . 

m=  1 


After  placing  these  equations  in  matrix  form  and  using  the  L-expanded  versions  of  u  and  v, 


x((k+l)L)  =  A1*  x(kL)  +  [  A^'^B  A'^2)B  ..  AB  B  ]  U(kL) 


(2.2) 


-  c  - 

D 

0 

0 

...  o  - 

CA 

CB 

D 

0 

...  0 

Y(kL)  = 

CA2 

x(kL)  + 

CAB 

CB 

D 

...  0 

U(kL) 

.ca^L 

1 

n 

> 

W 

CB  D  _ 

These  equations  have  the  form  of  a  time-invariant  discrete-ume  system  with  period  L: 

x((k+l)L)  =  Ae  x(kL)  +  Be  U(kL) 

,  k  =  0,l . 

Y(kL)  =  Ce  x(kL)  +  De  U(kL) 

The  matrices  (Ae,Be,Ce,De)  will  be  referred  to  as  the  single-bundle  L-expanded 
representation  of  (2.1),  where  “single-bundle”  refers  to  the  single-bundle  nature  of  the  inputs 
and  outputs.  Since  in  this  form  the  system  operates  on  blocks  of  input  values  and  produces  blocks 
of  output  values,  the  concepts  employed  above  are  often  called  block  processing  in  the  literature: 
see  Meyer  and  Burrus  (1975)  and  Burrus  (1972). 

The  L-expanded  representation  of  a  single-rate  discrete-time  system  with  inputs  and  outputs 
that  are  vectors  of  bundles  may  be  derived  using  manipulations  performed  in  finding  single-bundle 
L-expanded  representations.  Later  sections  will  show  that  multirate  discrete-time  systems  and 
multirate  sampled  continuous- time  systems  can  be  expressed  as  single-rate  discrete-time  systems, 
so  the  qualifier  “single-rate”  above  does  not  restrict  the  applicability  of  the  results  in  this  section. 

Let  u  and  y  be  vectors 


-ui  - 

■yi ' 

u  = 

U2 

and  y  = 

yi 

_um_ 

-yq- 

where  each  Uj,  j  =  l,...,m  and  each  y,,  i  =  l,...,q  are  bundles  of  signals,  and  all  bundles  are  at  the 
same  rate,  1/Tj.  Consider  a  single-rate  discrete-time  system  at  rate  \fT} : 

x((m+l)Tj)  =  A  x(mTp  +  B  u(mTj) 

,  m  =  0,1,...  .  (2.3) 

y(mTj)  =  C  x(mTj)  +  D  u(mTj) 

Let  L/Tj  €  N  and  B,  C,  and  D  be  partitioned  to  conform  with  the  bundles  in  u  and  y: 


"Cl" 

"  du  . 

dim 

B  =  [  bi  b2  ...  bm  ],  c  = 

: 

,  and  D  = 

-Cq- 

_  dqi 

■  dqm  _ 

II 


For  each  j  6  and  i  e  {l,...,q},  let  (Ae,bje,Cie,dije)  he  the  single-bundle  L-expanded 

representation  of  (A,bj,Ci,dij).  By  superposition, 

x((k+l)L)  =  Ae  x(kL)  +  Be  U(kL) 

Y(kL)  =  Ce  x(kL)  +  De  U(kL), 

where 

Cie  ^lle  —  *^lme 

Be  =  [  h)e  •  ■ '  bmg  ],  Cg  ,  and  Dg =  •  •  , 

_  Cqe  _  _  dqie  ...  dqme  _ 

and  Y(kL),  U(kL)  are  the  L-expanded  versions  of  y  and  u.  In  other  words,  (Ae,Be,Ce,De)  is  the 
L-expanded  representation  of  (2.3).  Be,  Ce,  and  De  are  partitioned  conformal  with  the 
expanded  bundles  composing  U(kL)  and  Y(kL).  As  in  the  case  of  signals,  L-expanded 
representations  may  be  expanded  again  by  thinking  of  (Ae,Be,Ce,De)  as  a  realization  for  a  single¬ 
rate  discrete-time  system  with  rate  1/L. 

For  a  system  where  the  inputs  and  outputs  have  multiple  bundles.  Finding  (Ae,Bc,Ce,De) 
by  merely  using  the  same  formulas  as  for  the  single-bundle  expansion  would  be  much  easier,  but 
this  results  in  values  over  time  L  of  each  signal  bundle  being  widely  scattered  throughout  U(kL) 
and  Y(kL).  By  performing  the  expansion  as  shown,  values  over  time  L  of  each  particular  signal 
bundle  are  adjacent  in  U(kL)  and  Y(kL).  This  greatly  facilitates  the  connection  of  expanded 
systems.  The  partition  of  U(kL)  and  Y(kL)  into  expanded  bundles  in  this  form  of  system 
expansion  reflects  the  designation  of  bundles  as  input  and  output  lines  on  a  block  diagram. 

The  expanded  representation  possesses  many  of  the  characteristics  of  the  original 
realization.  It  is  easily  verified  that  the  number  of  states,  stability,  reachability,  and  observability 
of  the  original  realization  are  all  preserved  by  expansion.  A  feature  of  expanded  representations 
that  complicates  certain  control  applications  is  that  for  a  system  with  no  direct  feedthrough. 
iA.B.C.O),  the  expanded  representation  will  have  De  nonzero  in  general.  Computationally, 
expansion  involves  only  multiplication,  composition,  and  storage  of  matrices.  For  large  1„  storing 
the  distended  matrices  of  the  expanded  representation  may  present  difficulties:  Bc  and  C0  each 
grow  linearly  with  1„  and  the  size  of  De  increases  as  the  square  of  1,. 

The  expanded  representation  of  a  subsystem  represents  a  step  toward  the  goal  of  finding  a 
time-invariant  description  of  a  multirate  system.  To  complete  the  process  of  finding  a  time- 
invariant  description  of  a  multirate  system,  the  expanded  representation  of  multirate  sampled 
continuous-time  systems  will  be  found,  and  methods  for  combining  expanded  representations  will 
be  developed. 


2.4  Expanding  Multirate  Sampled  Continuous-Time  Systems 


To  find  the  expanded  representation  of  a  multirate  sampled  continuous-time  system,  it  must 
be  discretized  first.  To  perform  discretization,  each  input  or  output  of  the  system  must  be  sampled 
at  only  one  rate.  Inserting  additional  samplers  that  do  not  alter  the  behavior  of  the  system  at 
strategic  locations  in  the  block  diagram  of  the  multirate  system  facilitates  discretization  of  the 
continuous-time  subsystems  and  will  also  simplify  the  operations  to  be  performed  in  Section  2.5. 
Before  discretizing  a  continuous-time  system,  it  is  assumed  that  the  following  procedure  is 
followed: 

Procedure  2.2: 

Assume  that  the  multirate  system  in  question  satisfies  Assumption  2.1  and  that  Procedure 
2.1  has  been  performed  The  subsystems  referred  to  in  the  following  steps  are  the  subsystems 
which  result  from  the  partitioning  required  by  Assumption  2.1(e). 

a.  Insert  a  sample  and  hold  on  each  input  and  output  of  every  discrete-time 
subsystem  in  the  block  diagram  (points  b,d,e,g,  and  i  in  Figure  2.3):  its  rate  is  that 
of  the  respective  output  or  input  of  that  subsystem. 

b.  Refer  to  both  an  output  of  a  continuous-time  subsystem  and  an  input  to 
the  entire  multirate  system  as  a  (system)  continuous  output  point.  For  each 
continuous  output  point  in  the  multirate  system  that  is  not  immediately  followed  by 

a  sampler  (points  a  and  c  in  Figure  2.3,)  let  {Tj . Tj )  be  the  set  of  periods  of  the 

first  sampler  encountered  on  each  path  leaving  that  continuous  output  point.  Insert 
into  the  block  diagram  a  sample  and  hold  with  period 

Tj  =  t(GCD  { p! . Pj })  =  pjt 

immediately  after  that  continuous  output  point 

c.  Refer  to  both  an  input  to  a  continuous-time  subsystem  and  an  output  of 
the  entire  multirate  system  as  a  (system)  continuous  input  point.  For  each 
continuous  input  point  in  the  multirate  system  that  is  not  immediately  preceded  by  a 
sampler  (points  h  and  j  in  Figure  2.3,)  let  {Ti,...,Tj)  be  the  set  of  periods  of  the 
last  sampler  encountered  on  each  path  leading  to  that  continuous  input  point.  Insen 
into  the  block  diagram  a  sample  and  hold  with  period 

Tj  =  t(GCD  { pi . pj))  =  pjt 

immediately  before  that  continuous  input  point. 

d.  Refer  to  an  output  of  a  summer  not  in  one  of  the  subsystems  composing 
the  multirate  system  (points  f,  g,  and  j  in  Figure  2.3)  as  a  (connection) 
continuous  output  point.  Repeat  the  actions  performed  for  system  continuous 


output  points  in  step  (b)  for  each  connection  continuous  output  point  in  the  multirate 
system. 

e.  Refer  to  an  input  of  a  branch  point  not  in  one  of  the  subsystems 
composing  the  multirate  system  (points  a,  e,  and  f  in  Figure  2.3)  as  a 
(connection)  continuous  input  point.  Repeat  the  actions  performed  for 
system  continuous  input  points  in  step  (c)  for  each  connection  continuous  input 
point  in  the  multirate  system. 

♦  ♦  ♦ 

The  satisfaction  of  parts  (e)  and  (0  of  Assumption  2.1  assures  the  existence  of  the  rates 
required  in  parts  (b)  and  (c)  of  Procedure  2.2.  All  possible  sources/destinations  of  signals  in  the 
multirate  system  are  thus  followed/preceded  by  samplers;  consequently,  the  rates  required  in  parts 
(d)  and  (e)  of  Procedure  2.2  will  exist.  The  application  of  Procedure  2.2  does  not  alter  the 
behavior  of  the  system;  each  added  sampler  updates  at  least  at  the  times  that  the  samplers  feeding 
from/to  it  update.  By  definition  of  Tj  (see  parts  (b)  and  (c)  of  Procedure  2.2,)  if  a  sampler  with 
period  T„  i  s  ( l,...,j),  samples,  the  sampler  with  period  Tj  also  samples;  t  =  nT„  ne?  implies 
t/Tj  =  npit/pjT  =  n(pi/pj)  e  Zf 


Figure  2.3.  C\  and  C2  are  continuous-time  subsystems,  and 
Di  and  D;  are  discrete-time  subsystems. 


Example  2.3; 

Let  Di  and  D2  in  Figure  2.3  be  single  rate  with  rates  1/T 1  and  I/T2, 
respectively.  The  application  of  Procedure  2.2  to  the  multirate  system  in  Figure  2.3 
proceeds  as  follows: 

a:  Insert  samplers  of  period  Tb  =  T<j  =  Te  =  Tj  at  points  b.d.  and  e  and 
samplers  of  period  Ts  =  T,  =  T2  at  points  g  and  i. 
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b:  Insen  a  sampler  of  period  Ta  =  x(GCD  (p4,pb))  at  point  a  and  a  sampler 
of  period  Tc  =  x(GCD  {p5,pg})  at  point  c. 

c:  Insen  a  sampler  of  period  Tj  =  x(GCD  { P6>Pi ) )  at  point  j  and  a  sampler 
of  period  Th  =  Te  at  point  h. 

d:  Insen  a  sampler  of  period  Tf  =  x(GCD  { Ps.Pg} )  at  point  f.  Points  g  and 
j  are  also  connection  continuous  output  points,  but  samplers  are  already  present  at 
these  points. 

e:  Points  a,e,  and  f  are  connection  continuous  input  points,  but  samplers 
are  already  present  at  these  points. 


After  applying  Procedure  2.2,  each  input  and  output  line  of  a  continuous-time  subsystem 
on  the  block  diagram  should  be  designated  as  a  signal  bundle.  The  last  sampler  an  input  bundle 
passes  through  or  the  first  sampler  an  output  bundle  passes  through  determines  the  rate  assigned  to 
that  bundle  during  the  discretization  process. 

Discretization  of  a  multirate  sampled  continuous-time  system  is  also  considered  in  Araki 
and  Yamamoto  (1986),  from  a  perspective  quite  different  from  the  one  taken  here;  Araki  and 
Yamamoto  (1986)  employ  expanded  states,  as  well  as  expanded  signals. 

Let  the  continuous- time  system 

x=Ax  +  Bu,  y  =  Cx  +  Du  (2.4) 

have  the  input 

~ui  “ 
u2 


where  each  u;  is  a  signal  bundle  sampled  and  held  at  rate  1/Tj,  i  =  l,...,m,  and  B  is  partitioned  to 
be  compatible  with  u:  B  =  [bi  t>2  ...  bnJ.  Let  y  be  partitioned  into  q  bundles,  with  each  bundle  Vj 
sampled  with  period  Tj,  j  =  l,...,q;  let  C  and  D  be  partitioned  accordingly: 


,  C  = 


,  and  D  = 


dql  . . .  dqm 


To  discretize  this  system,  let  L  =  LCP  {Tj,...  *Tm>Ti ,...  ,Tq ) ,  x(kL)  €  Rn  denote  the  state  of  (2.4) 
at  time  kL,  and  Y(kL)  and  U(kL)  denote  the  L-expanded  versions  of  v  and  u.  Let 

Si  =  T/T„  sj  =  T/Tj,  <D(w)  =  eAw,  and  T(w)  =  J*  eA(w0dt  . 


S3 

ft 
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Then  a  straightforward  but  tedious  calculation  gives 


x((k+l)L)  =  E  x(kL)  +  F  U(kL) 
Y(kL)  =  G  x(kL)  +  H  U(kL), 


where 


1.  E  =  d>(L). 


2.  F  =  [fi  f2 ...  fnJ,  partitioned  compatible  with  U(kL),  where 

fi  =  [OCriCsi-l))^^^  I ...  1  <D(Ti)r(Ti)bi  I  r(Ti)bi]  ,  i  =  l,...,m. 


g2  CjO(Tj) 

G  =  is  partitioned  as  Y(kL),  where  gj  = 


CjdKT/Sj-l)) 


,  j  =  1»— »q. 


hi i  ...  hj 


L  hql  • • •  hqm 


is  partitioned  according  to  Y(kL)  and  U(kL),  and  each  subblock  h;j,  i  =  l,...,q,  j  =  l,...,m,  is 
composed  of  blocks  (h;j)rv,  r  =  l,...,Si,  v  =  l,...,Sj,  each  having  the  same  dimension  as  d;j  and 
where 


0  ,  if  (r-l)Ti  <  (v-l)Tj 


(hIJ)rv=<  dij  +  c;r((r-l)Ti  -  (v-l)Tj)bj  ,  if  (v-l)Tj  <  (r-l)Ti  <  vT;  . 


CiO((r-l)Ti  -  vTprCTjlbj  ,  if  vTj  <  (r-l)T, 


From  (2.5)  and  points  1  through  4  above,  the  discretized  multirate  sampled  continuous-time 
subsystem  is  a  time-invariant  single-rate  discrete-time  system  with  inputs  and  outputs  that  are 
vectors  of  L-expanded  bundles.  Thus,  the  methods  in  Section  2.3  apply  to  further  expansion  of 
the  discretized  continuous- time  system. 

Despite  the  formidable  appearance  of  the  discretized  equations,  the  required  calculations  can 
be  performed  using  any  software  package  capable  of  discretizing  a  continuous-time  state  equation 
at  a  single  rate.  The  quantity  T(w)B  is  the  “B”  matrix  resulting  from  discretizing  (A,B)  at  rate  !/'w. 
The  required  arguments  of  d>(.)  and  f(.)  are  all  integer  multiples  of  x.  For  some  values  of  i,  j,  r. 


V  V-V.  v 


and  v,  the  arguments  of  <t>(.)  and  T(.)  will  coincide,  so  not  every  occurrence  of  OC.)  and  T(.)  in  the 
above  equations  entails  an  additional  calculation. 

A  modification  to  Procedure  2.2,  which  greatly  simplifies  the  discretization  of  multirate 
sampled  connnuous-time  subsystems,  is  to  place  sample/holds  of  period  x  on  the  block  diagram  at 
the  inputs  and  outputs  of  each  continuous-time  subsystem,  leaving  Procedure  2.2  otherwise 
unchanged.  These  added  samplers  leave  the  behavior  of  the  system  unaltered  for  the  same  reason 
as  that  given  following  Procedure  2.2.  This  modification  has  the  advantage  that  discretization  can 
be  earned  out  uniformly  at  rate  t,  but  it  may  result  in  inputs  and  outputs  of  unacceptably  large 
dimensions  after  expansion 


Example  2  4: 


Let  the  continuous-time  subsystem  in  Figure  2.4  have  the  state  equation 

x(t)  =  A  x(t)  +  B  u(t) 

ryiwi  rcil  r^n  .. 


yi(t)i  fen  rdf 

=  x(t)  +  u(t). 

y2(t)J  Lc2J  Ld2_ 


Additional  samplers  do  not  need  to  be  added  to  discretize  this  system.  Following 
the  steps  above,  let 

O(w)  =  eAw,  and  T(w)  =  JQW  eA(w-t)dt . 


x((n+l)6t)  =  <P(6x)  x(n(6t))  +  [<D(3x)r(3x)B  R3x)B] 


u(n(6t)) 

u(n(6x)+3x) 


yiCn(6x)) 

yi(n(6x)+x) 

yi(n(6x}+2x) 

y!(n(6x)+3t) 

yi(n(6t)+4i) 

yi(n(6x)+5x) 

v2(n(6x)) 

y2(n(6x)+2x) 

,y2(n(6x)+4x). 


ciO(x) 

ciO(2x) 

Ci<D(3x) 

CiO(4x)  x(n(6x)) 

ciO(5x) 

C2 

C2<1>(2x) 

.c20(4x)_ 


di+ciT(x)B  0 

di+ciT(2x)B  0 

CiT(3x)B  dj 

Ci<E>(x)T(3x)B  di+CiRxlB 

CiO(2x)r(3x)B  di+c,r(2x)B 
d2  0 

d2+c2r(2x)B  0 

c20(x)F(3x)B  d2+c2r(x)B 


u(n(6x)) 

u(n(6x)-f3r) 


Note  that  all  of  the  quantities  above  could  be  obtained  by  normal  discretization  with 
periods  x,  2x,  and  3x.  Also,  if  the  modification  to  Procedure  2.2  noted  above  were 


used,  the  resulting  discretized  system  would  have  6  inputs  and  12  outputs. 

♦  ♦  ♦ 


Figure  2.4.  A  multirate  sampled  continuous-time  subsystem. 

2.5  Connection  Matrices 

The  concept  of  a  connection  matrix  facilitates  the  interconnection  of  subsystems  whose 
inputs  and  outputs  are  expanded  bundles  of  signals  which  pass  through  summers,  samplers,  and 
branch  points  between  the  subsystems.  An  individual  connection  matrix  serves  to  connect  a  signal 
bundle  at  one  rate  to  a  signal  bundle  at  another  rate  after  expansion.  The  form  of  the  connection 
matrix  depends  only  on  the  rates  of  the  bundles  being  connected  and  the  amount  of  time  over 
which  the  signals  were  expanded. 


Figure  2.5.  A  direct  connection  between  subsystems. 

Consider  a  single-rate  subsystem,  £i,  with  a  bundle  of  r  outputs  connected  directly  to  a 
single-rate  subsystem  with  a  bundle  of  r  inputs,  £2-  as  shown  in  Figure  2.5.  Let  L/Tj  =  lj, 
L/T;  =  1;,  and  lj,  I2  €  N.  Since  h  *  I2  in  general,  the  number  of  outputs  of  the  L-expanded 
description  of  £1  will  not  equal  the  number  of  inputs  of  the  L-expanded  description  of  £2-  The 
mechanism  which  serves  to  connect  Yi(kL)  to  U2(kL)  is  the  sampling  £2  performs  on  its  input: 

U2<kL+(i-l)T2)  =  yi(kL+(j-l)Ti)  if  (j-l)Ti  <  (i-l)T2  <  jTt  , 

where  i  €  { 1  ,...,12 }  and  j  €  1 1  ,...,lj }.  This  relationship  can  be  expressed  concisely  by  defining 
Q(L.r,l:2),  the  connection  matrix  for  L-expanded  bundles  of  r  signals  from  period  Tj  to  period 
T:,  such  that  LT2(kL)  =  Q(L,r,l:2)Yi(kL).  Q(L,r,l:2)  is  an  rl2  x  rlj  matrix  of  L  x  lj  rxr  blocks, 
w  here  the  i,j&  block  is  given  by 

r  I-  ,  if  (j- 1  )pi  <  (i- 1  )p2  <  jp  1 


100  000 
100  000 
010  000 
010  000 
001  000 
000  100 
000  100 
000  010 
000  010 
000  001 

♦  ♦  ♦ 

A  fact  which  is  easily  proven  by  induction  is  that  if  M/T]  e  N  and  M/T2  e  N,  then  for 
some  me  N,  M  =  mt(LCM  { Pi  ,p2 } )-  In  words,  any  M  admissible  for  expanding  signals  with 
periods  Ti  and  T2  is  an  integer  multiple  of  the  least  common  period  of  Ti  and  T;.  Let 
T(l,2)  =  t(LCM  {pt,p2}).  Since  U2(nT(  1,2))  =  Q(T(l,2),r,l:2)Yl(nT(l,2)),  this  connection 
matrix  may  be  thought  of  as  a  discrete-time  system  with  rate  1/T(  1 ,2),  input  Yi(nT(l,2)),  output 
U2(nT(l,2)),  and  representation  (0,0,0,Q(T(l,2),r,l:2)).  After  finding  the  M-expanded 
representation  for  this  system, 

U2(kM)  =  diag[Q(T(  1 ,2),r,  1 :2),Q(T(  1 ,2),r,  1 :2) . Q(T(l,2),r,l:2)]  Y^kM): 

Q(T(l,2),r,l:2)  appears  m  times  along  the  diagonal.  Thus, 

Q(M,r,  1 :2)  =  diag[Q(T(  1 ,2),r,  1 :2),Q(T(  1 ,2),r,  1 :2),...,Q(T(  1 ,2),r,  1:2)1. 

This  property  greatly  reduces  the  effort  required  to  compute  and  store  connection  matrices  for  large 
m.  In  Example  2.5,  m  equals  2. 

Connection  matrices  can  be  systematically  placed  so  that  a  collection  of  L-expanded 
subsystems  is  connected  according  to  a  given  block  diagram.  Since  separate  inputs  and  outputs 
drawn  on  a  block  diagram  were  defined  as  distinct  bundles  of  signals,  connection  matrices  can 
effect  direct  connections  between  subsystems.  Under  the  assumption  that  Procedure  2.2  has  been 
applied  to  the  block  diagram,  rules  for  placing  connection  matrices  when  the  bundles  of  signals 
pass  through  samplers,  summers,  and  branch  points  between  the  subsystems  are  easily  stated. 

Define  a  link  as  any  path  from  one  sampler  to  another  on  the  block  diagram  which  may 
pass  through  summers  and  branch  points  but  not  through  subsystems  or  other  samplers.  As  a 
result  of  applying  Procedure  2.2,  in  the  region  between  the  subsystems  a  sampler  follows  each 
summer  and  precedes  each  branch  point.  Thus,  a  link  passes  through  at  most  one  summer  and  one 


branch  point.  The  steps  necessary  to  connect  a  collection  of  L-expanded  subsystems  together 
according  to  a  given  block  diagram  are  detailed  below. 

Procedure  2.3: 

Assume  that  all  subsystems  and  signals  in  the  block  diagram  are  in  L-expanded  form  and 
that  Procedure  2.2  has  been  applied  to  the  block  diagram.  For  each  link  in  the  system  to  be 
connected: 

a.  Find  Q(L,r,f:t),  where  the  r  signals  on  the  link  flow  from  a  sampler  with 
period  Tf  to  a  sampler  with  period  Tt. 

b.  Insert  the  connection  matrix  in  the  link  according  to  the  type  of  link  as 
follows: 

i.  If  the  link  traverses  a  branch  point  and  a  summer,  place 
the  connection  matrix  after  the  branch  point  and  before  the  summer. 

ii.  If  the  link  crosses  only  a  branch  point,  place  the 
connection  matrix  after  the  branch  point 

iii.  If  the  link  traverses  only  a  summer,  insert  the  connection 
matrix  before  the  summer. 

iv.  If  none  of  the  above  apply,  the  link  crosses  no  summers 
or  branch  points,  and  the  connection  matrix  can  be  placed  anywhere 
on  the  link. 

♦  ♦  ♦ 

As  an  aid  to  justifying  the  connection  matrix  locations  specified  in  Procedure  2.3,  note  that  their 
placement  ensures  that  the  dimensions  of  all  inputs  of  each  summer  and  of  all  outputs  of  each 
branch  point  in  the  block  diagram  of  the  expanded  system  are  identical. 


I 
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Example  2.6: 

Figure  2.6  shows  the  placement  of  connection  matrices  in  the  multirate 
system  in  Figure  2.3  after  Procedure  2.2  has  been  performed.  The  numbers  1,  2, 
3,  and  4  in  Figure  2.6  denote  a  link  of  type  i,  ii,  iii,  or  iv,  respectively,  as  described 
in  Procedure  2.3(b). 


2.6  The  T-expanded  Representation  of  a  Multirate  System 


The  preceding  sections  collectively  provide  a  means  of  obtaining  a  time-invariant  input- 
output  or  state  space  description  of  a  multirate  system,  which  will  be  referred  to  as  the 
T-expanded  representation.  The  reader  may  wish  to  keep  the  following  general  procedure  in 
mind  as  a  summary  of  the  previous  sections. 

Procedure  2.4: 

a.  Starting  with  a  block  diagram  of  the  system,  verify  that  Assumption  2. 1 
is  satisfied.  The  partitioning  required  in  Assumption  2.1(e)  may  entail  a  trial  and 
error  approach.  For  the  remaining  steps  of  this  procedure,  “subsystem”  refers 
specifically  to  the  subsystems  obtained  as  a  result  of  this  partitioning. 

b.  Assign  labels  to  all  sampling  periods  in  the  system,  and  apply  Procedure 
2.1  to  find  T. 

c.  Find  a  state  space  realization  for  each  continuous-time  and  discrete-time 
subsystem. 

d.  Apply  Procedure  2.2  to  the  block  diagram.  Procedure  2.2(b)  may  be 
performed  in  parallel  with  Procedure  2.2(c),  and  Procedure  2.2(d)  and  Procedure 
2.2(e)  may  be  performed  in  parallel.  However,  do  not  reverse  the  order  of  these 
two  pairs  of  steps. 

e.  Identify  the  rates  at  which  the  inputs  and  outputs  of  continuous-time 
subsystems  are  sampled  and  discretize  these  subsystems  as  described  in  Section 
2.4. 

f.  Designate  lines  on  the  block  diagram  representing  inputs  and  outputs  of 
subsystems  and  inputs  and  outputs  of  the  entire  multirate  system  as  bundles  of 
signals.  Partition  the  state  space  realization  of  each  subsystem  accordingly. 

g.  Find  the  T-expanded  representation  of  each  subsystem  as  detailed  in 
Section  2.3. 
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h.  Apply  Procedure  2.3  and  then  delete  all  of  the  samplers  from  the  block 
diagram. 

i.  Note  that  the  block  diagram  obtained  in  step  (h)  is  a  block  diagram  of 
interconnected  time-invariant  discrete-time  systems  with  period  T.  The  states  at 
times  kT  of  the  subsystems  in  this  block  diagram  and  the  input-output  behavior  of 
the  entire  block  diagram  as  a  system  are  identical  to  those  of  the  original  multirate 
system  in  step  (c)  (after  decomposing  the  expanded  signals  into  their  components.) 
Standard  methods  may  now  be  applied  to  reduce  this  block  diagram  to  a  single  state 
space  equation  or  transfer  function  matrix. 

♦  ♦  ♦ 

Completion  of  Procedure  2.4(i)  will  in  general  involve  a  matrix  inversion  to  find  the  state 
space  description  of  a  feedback  structure.  Under  certain  conditions,  this  inverse  may  not  exist; 
hence,  a  state  space  representation  for  such  a  feedback  structure  cannot  be  found.  As  an  example, 
the  seemingly  innocuous  system  in  Figure  2.7  composed  of  two  samplers  and  two  continuous- time 
unity  gain  blocks  has  no  state  space  description  for  any  choices  of  Ti  and  T2.  This  difficulty  is 
intimately  related  to  the  well-posedness  problem  encountered  in  transfer  function  descriptions  of 
composite  systems.  See  Chen  (1984)  for  a  discussion  of  this  problem.  If  the  multirate  system  as 
modeled  by  the  block  diagram  in  Procedure  2.4(a)  has  the  property  that  each  closed  path  on  the 
block  diagram  passes  through  at  least  one  subsystem  which  has  no  direct  feedthrough,  then  it  is 
suspected  that  the  difficulties  described  above  will  not  be  encountered. 


Figure  2.7.  An  ill-posed  multirate  system. 

2.7  Inclusion  of  Periodic  Subsystems 

A  slight  extension  of  Procedure  2.4  permits  the  inclusion  of  single-rate  periodic  discrete¬ 
time  subsystems  in  the  multirate  system.  Consider  a  single-rate  discrete-time  subsystem  with 
single-bundle  inputs  and  outputs  and  a  realization  which  is  time  varying,  but  periodic.  For  all 
n  €  Z+,  let 

KOi+DTO  =  I(nT,)  =  (A(nTi),B(nTi),C(nT1),D(nT1)). 

Thus,  I(nTi)  is  periodic  with  period  rTj.  As  demonstrated  in  Meyer  and  Burrus  (1975),  repeated 


iteration  from  an  initial  time  of  t  =  0  of  the  state  equation  governing  I(nTi), 
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x((n+l)Ti)  =  A(nTj)  x(nTi)  +  B(nTi)  u(nTi) 
y(nTi)  =  CfnTO  x(nTi)  +  D(nTi)  u(nT0, 


yields  a  time-invariant  state  equation  with  rate  1/rTi: 


x((k+l)rTO  =  Ae  x(k(rTi))  +  Be  U(k(rT0) 
YfkfrTO)  =  CexfkdTO)  +  DeUfkfrTO). 


Denoting  A(nTi),  B(nTi),  etc.  as  A(n),  B(n),  etc.  for  brevity, 


Ae  =  A(r-l)A(r-2)...A(l)A(0), 


Be=  [A(r-1)...A(1)B(0)  I...I  A(r-l)A(r-2)B(r-3)  I  A(r-l)B(r-2)  I  B(r-l)], 


C(0) 

C(1)A(0) 


C,  =  C(2)A(1)A(0) 


.C(r-l)A(r-2)...A(0) 


and  De  is  composed  of  r  x  r  of  the  blocks  djj,  where 

"  0  ,  if  i  <  j 
=  <  D(i-l)  ,  if  i  =  j 
ij  |  C(i-l)B(j-l)  ,  ifi=j+l 

.  C(i-l)A(i-2)...A(j)B(j- 1)  ,  if  i  >  j+1 


Thus,  (Ae,Be,Ce,De)  serves  as  a  single-bundle  rTi  -expanded  representation  of  I(nTi).  This  time- 
invariant  representation  can  be  further  expanded  using  the  technique  in  Section  2.3.  If  Procedure 
2.4  were  applied  to  a  multirate  system  containing  this  periodic  subsystem,  the  only  alteration  of 
Procedure  2.4  required  would  be  to  include  both  Ti  and  rTj  in  Procedure  2.4(b)  to  find  T.  In 
analogy  with  Section  2.3,  the  expanded  representation  of  a  periodic  discrete-time  subsystem 
S(nTi)  with  inputs  u(nTi)  and  outputs  v(nTi)  that  are  vectors  of  bundles  can  be  expressed  in 
terms  of  single-bundle  expanded  representations  by  first  partitioning  B(nTi),  £(nT0,  and  QinTi ; 
to  conform  with  the  bundles  composing  u(nTi)  and  y(nTi). 

A  multirate  periodic  discrete-time  subsystem  may  be  specified  in  Tj-expanded  form 
by  both  specifying  A(kTi),  B(kTi),  QkTi),  and  D(kTi)  in  its  state  equation. 


x((k+l)Tj)  =  A(kTj)  x(kTi)  +  B(kTO  U(kT,) 
Y(kTj)  =  C(kT0  x(kTi)  +  D(kTi)  U(kT0, 


where  Y(kTi)  and  U(kTi)  are  Tj-expanded  versions  of  vectors  of  bundles  of  signals  and  A(kTj), 
B(kTt),  C(kTi),  and  D(kTi)  are  rTj -periodic,  and  by  specifying  the  partition  of  Y(kTj)  and 
U(kTi)  into  expanded  bundles  and  the  rates  associated  with  each  expanded  bundle.  Although 
discretization  of  a  rTi-periodic  continuous-time  system  with  inputs  and  outputs  sampled  at  multiple 
rates  would  yield  a  state  equation  such  as  (2.7),  the  state  equation  (2.7)  may  not  be  the  result  of 
discretizing  or  expanding  any  subsystem.  In  the  case  r  =  1,  (2.7)  may  represent  the  behavior  of  a 
computer  program  performing  concurrent  tasks  at  multiple  rates.  In  Chapter  5,  the  controllers 
designed  will  take  the  form  of  (2.7).  The  subsystem  (2.7)  may  be  expanded  over  its  period,  rTi, 
and  then  further  expanded  over  time  T  during  the  execution  of  Procedure  2.4. 

The  point  of  the  extensions  of  Procedure  2.4  given  above  is  not  to  belabor  specific 
examples  but  to  indicate  a  general  property  of  Procedure  2.4.  If,  by  some  means,  a  time- invariant 
discrete-time  description  of  a  subsystem  in  terms  of  expanded  inputs  and  outputs  can  be  obtained, 
Procedure  2.4  applies  to  a  multirate  system  in  which  this  subsystem  appears  as  a  component.  In 
this  instance.  Procedure  2.4  must  only  be  modified  by  including  the  period  of  the  time-invariant 
description  of  this  subsystem  in  the  set  of  sampling  periods  considered  when  finding  the  least 
common  period  of  the  entire  multirate  system,  T. 

Example  2.7: 

Suppose  that  a  T-expanded  representation  is  desired  for  the  multirate  system 
shown  in  Figure  2.8,  where  Ti  =  1  sec,  T2  =  0.5  sec,  and  g(nTj)  is  a  discrete-time 
periodic  gain  with  period  T3  =  2Ti: 

f  gi  ,  if  n  is  even 

g(nTl)=  .f  .  • 

I  g2  ,  if  n  is  odd 


Figure  2.8.  A  multirate  system. 


Applying  Procedure  2.4: 

a:  Considering  the  double  integrator  and  the  discrete-time  gain  as 
subsystems,  parts  (a),  (e),  and  (f)  of  Assumption  2.1  are  satisfied.  It  will  be 
assumed  that  the  rest  of  Assumption  2.1  is  satisfied  as  well  (with  the  exception  that 
a  periodic  discrete-time  subsystem  is  present.) 
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b:  Application  of  Procedure  2.1  yields  x  =  0.5  sec,  pi  =  2,  p2  =  1,  P3  =  4, 
and  T  =  2  sec. 

c:  A  state  space  realization  for  the  periodic  gain  is  (0,0,0,g(nTj)),  and  a 
state  space  realization  for  the  double  integrator  is 

(  T  0  1  1  r  0  1  A 


f  0  11  r  01 

,  ,  [  1  0],  0 
L  o  o  J  L  i  J 


d:  If  the  modified  version  of  Procedure  2.2  (wherein  continuous-time 
subsystems  are  sampled  with  period  t  (=  T2))  is  applied  to  the  system,  Figure  2.9 
results. 


■  m 

u  i,  y-' — ► 

T.  V'  T, 

(1/s) 

- - 

T, 

1  +  A  2 

2 

m 

1 - N_ 

g(nT, ) 

◄ — nJ 

(2.6), 


Figure  2.9.  The  system  in  Figure  2.8  with  inserted  samplers, 
e:  Discretizing  the  double  integrator  with  period  T2  gives 

r  1  1/2  1  n/81 

Ad  —  ,  Bd  =  ,  Cd  =  [  1  0  ]  ,  and  Dd  =  0. 

L  0  1  J  L1/2J 

f:  This  part  is  trivial  since  ail  inputs  and  outputs  are  drawn  as  single  lines, 
g:  Denote  the  T-expanded  representation  of  the  periodic  gain  by  K.  From 

K=[8'  °\ 

.  0  g2  . 


The  T-expanded  representation  of  (Ad,Bd,Cd,Dd)  can  be  found  using  (2.2)  with 
It  =  T/T  2  —  4: 

r  1  2  1  f  7/8  5/8  3/8  1/8  1 
Ae=  I.  Be-  _  I. 


r  1  2 1 

r  7/8 

5/8 

3/8 

1/8 1 

,  B,= 

L  0  1 J 

L 1/2 

1/2 

1/2 

1/2 . 

1 

0  -j 

-  0 

0 

0 

0 

1 

1/2 

1/8 

0 

0 

0 

, and  De  = 

1 

1 

3/8 

1/8 

0 

0 

1 

3/2  J 

5/8 

3/8 

1/8 

0 

h:  LetQj  =Q(T,1,1:2),  Q2  =  Q(T,  1,2:1).  and  Q3  =  Q(T.  1.2:2).  Then 

25 


-  1  0  0  0  - 

-  Q:  =  0  0  10  ’  anc*  Q3  =  k- 


After  inserting  connection  matrices  and  removing  samplers,  the  block  diagram  of 
the  T-expanded  representation  of  the  system  appears  as  in  Figure  2.10. 


U(mT) 


Y(mT) 


(A  ,B  ,C  ,D  1 

e  e  e  e' 


Figure  2.10.  Block  diagram  of  the  T-expanded  representation. 

i:  The  system  in  Figure  2.10  can  be  reduced  to  a  single  state  space 
equation,  the  T-expanded  representation  of  the  system  in  Figure  2.8: 


"  4+6gi+gig2+2g2  4+g;  -i  r  6+g;  1  " 

4  2  4  2  f  u(mT) 

x((m+l)T)  =  x(mT)  + 

gtg2+2gi+2g2  2+g;  ,  Lu(mT+T]), 


y(mT)  - 
y(mT+T2) 
y(mT+2T2) 
y(mT+3T2). 


10-1  r  0  0 

1  1/2  1/8  0 

x(mT)  + 

11  1/20 


1  3/2 


u(mT)  1 
uin.T+Ti)J 


1  1/8 


2.8  Modular  Expansion  of  Multirate  Systems 

A  fundamental  difficulty  encountered  in  executing  Procedure  2.4  is  that  the  T-expanded 
representation  of  the  subsystems  may  involve  matrices  of  large  dimensions.  A  simple  example 
demonstrates  this  problem  and  suggests  a  method  of  circumventing  it.  referred  to  as  modular 
expansion. 

Consider  the  multirate  system  with  SISO  time-invariant  subsystems  in  Figure  2.1 1.  where 
T[  =  2,  T;  =  1.  T3  =  25,  and  T  =  50.  Straightforward  application  of  Procedure  2.4  yields  .1 


iw'i  ,  **»  ‘V*1.  * , 


T-expanded  representation  for  £1  with  25  inputs  and  outputs  and  a  T-expanded  representation  for 
£2  with  50  inputs  and  outputs.  Reduction  of  the  feedback  connection  of  the  expanded 
representations  to  a  single  state  equation  involves  the  multiplication  and  addition  of  large  matrices 
and,  more  importantly,  the  inversion  of  two  25  x  25  matrices. 


Figure  2.1 1.  A  simple  multirate  system. 


As  an  alternative  method,  consider  ignoring  £3  for  the  moment  and  finding  an  expanded 
representation  for  the  feedback  connection  of  Li  and  £2  over  the  least  common  period  ot  Tj  and 
T;  using  Procedure  2.4.  Let  Tf  =LCP  (Ti,T2)  =  2.  Calculation  of  the  Tf-expanded 
representation  for  the  feedback  connection  involves  much  smaller  matrices  and.  in  particular, 
requires  the  inversion  of  two  1  x  1  matrices.  This  savings  in  computation  will  be  worthwhile  if  the 
T-expanded  representation  of  £f  exists,  enabling  the  T-expanded  representations  of  £f  and  £3  to  be 
combined.  A  Tf-cxpanded  representation  can  be  T-expanded  if  T/Tf  e  N.  Since  Theorem  2.3 
asserts  that  T/Tf  =  LCP  { T i  T2,T3)/LCP  {Ti.TiJ  e  N,  the  T-expanded  representation  of  £t  can 
be  computed.  This  then  is  the  essence  of  modular  expansion. 

There  is  substantial  freedom  in  the  steps  taken  when  performing  modular  expansion  on  a 
multirate  system.  The  following  procedure  helps  explain  the  notion  of  modular  expansion  and 
should  not  be  interpreted  as  the  only  means  by  which  modular  expansion  can  be  performed. 
Procedure  2.5:  (a  modular  expansion  scheme) 

a.  Carry  out  steps  (a)  through  (f)  of  Procedure  2.4. 

b.  Focus  attention  on  a  collection  of  subsystems  by  drawing  a  path  on  the 
block  diagram  enclosing  one  or  more  subsystems  which  does  not  separate  any 
enclosed  continuous  output/input  point  from  the  sampler  immediately 
following/preceding  it  (review  Procedure  2.2  for  the  meaning  of  these  terms.  )  Call 
the  enclosed  portion  of  the  system  the  “current  collection.”  Set  i  =  0  and  P(0)  =  1 

Let  (T0.i,To.2 . T0.m(o)}  be  the  set  of  m(0)  sampling  periods  and  periodic 

discrete-time  subsystem  periods  in  the  current  collection. 


c.  Set  P(i+1)  =  LCM  (P(i),  pu . Pi,m(i)}  and  T(1+i)fo  =  P(i+l)x.  Find 

the  ^i+D.o-expanded  representation  of  each  subsystem  in  the  current  collection  and 
insert  connection  matrices.  By  standard  techniques,  find  a  single  discrete-time  state 
equation  with  rate  l/T^+i^o  to  describe  the  current  collection.  When  finding  this 
state  equation,  ignore  samplers  on  the  block  diagram  that  are  between  connection 
matrices,  subsystems,  summers,  or  branch  points  in  the  current  collection  but  retain 
samplers  on  the  periphery  of  the  current  collection. 

d.  Seti  =  i+1. 

e.  If  the  current  collection  is  the  entire  multirate  system,  stop.  Otherwise, 
draw  a  path  enclosing  the  current  collection  and  possibly  other  subsystems, 

following  the  same  restrictions  as  in  step  (b)  above.  Let  . Tiim(i)}  be 

the  m(i)  additional  sampling  periods  and  periodic  discrete-time  subsystem  periods 
enclosed  by  this  path.  Call  the  portion  of  the  system  enclosed  by  this  path  the 
current  collecdon  and  return  to  step  (c). 

♦  ♦  ♦ 

Procedure  2.5  proceeds  by  repeatedly  adding  subsystems  to  a  single  collection  of 
subsystems.  Modular  expansion  may  in  general  start  with  several  collections  scattered  throughout 
the  block  diagram  and  repeatedly  add  subsystems  and  collections  of  subsystems  to  these 
collections.  The  restrictions  in  Procedure  2.5(b)  on  the  items  included  in  a  collection  must  be 
obeyed  during  the  process  of  enlarging  a  collection.  Each  time  a  new  collection  is  formed,  that 
collection  is  expanded  and  then  reduced  to  a  single  state  equation  as  in  Procedure  2.5(c).  The 
process  terminates  when  a  single  collection  contains  the  entire  multirate  system.  The  principle  that 
allows  collections  of  subsystems  to  be  expanded  and  combined  in  such  a  variety  of  ways  is  that,  as 
a  simple  induction  argument  and  Theorem  2.3  show,  the  least  common  period  of  the  union  of  a 
collection  of  sets  of  sampling  periods  is  an  integer  multiple  of  the  least  common  period  of  each  set 
in  that  collection. 

A  trade-off  between  the  effort  required  to  combine  expanded  subsystems  and  the  effort 
expended  in  expanding  a  representation  corresponding  to  a  combination  of  subsystems  complicates 
the  question  of  which  modular  expansion  strategy  results  in  the  least  computational  effort  for  a 
given  multirate  system.  Two  general  heuristics  for  modular  expansion  can  be  offered.  To  avoid 
inverting  large  matrices,  expand  subsystems  in  a  feedback  connection  and  obtain  a  representation 
describing  this  feedback  configuration  before  including  it  in  a  larger  collection  of  subsystems.  By 
choosing  collections  of  subsystems  that  have  fewer  input  and  output  lines,  the  dimensions  of  the 
required  “B”,  “C”,  and  “D”  matrices  are  reduced. 


CHAPTER  3 


PERIODIC  REPRESENTATION  OF  MULTIRATE  SYSTEMS 


Representing  the  behavior  of  a  multirate  discrete-time  system  with  period  T  over  intervals 
of  time  of  length  M,  where  T/M  e  Q,  results  in  a  description  which  is  time  varying  and  periodic. 
Such  a  description  will  be  referred  to  as  the  M-varying  representation  of  the  system.  The 
motivation  for  examining  the  M-varying  representation  is  twofold.  Not  only  does  it  provide  a 
periodic  representation  of  the  multirate  system,  but  with  M  =  t,  the  M-varying  representation 
finds  use  theoretically  in  investigating  the  interperiod  behavior  of  a  multirate  system  by  examining 
its  T-expanded  representation.  Three  intrinsic  difficulties  detract  from  the  practical  usefulness  of 
the  M-varying  representation  presented  here.  In  general,  the  M-varying  representation  has  more 
states  than  the  time-invariant  description,  is  difficult  to  calculate,  and  involves  padding  with  false 
inputs  and  outputs  to  prevent  the  number  of  inputs  and  outputs  from  varying  with  time. 

Many  possibilities  exist  for  the  periodic  representation  of  a  multirate  system.  In  the 
M-varying  representation,  state  transitions  occur  at  regular  intervals.  A  representation  like  the  one 
presented  in  Kalman  and  Bertram  (1959)  in  which  state  transitions  occur  only  at  the  times  of 
certain  events,  such  as  a  sampling  event  or  a  state  transition  of  a  discrete-time  element,  is  perhaps 
more  appealing  from  a  practical  point  of  view  due  to  reduced  storage  requirements.  However,  the 
M-varying  representations  presented  here  are  valid  periodic  representations  of  the  multirate  system 
and  are  easily  described  in  detail. 


3.1  The  x-varying  Representation 

In  Chapter  2,  the  system  was  represented  over  an  interval  of  time  which  was  an  integer 
multiple  of  the  period  of  each  sampler.  Inputs  and  outputs  were  implicitly  “held”  by  retaining  their 
values  as  components  of  time  expanded  signals.  The  t-varying  representation  of  a  system  does  not 
permit  this  luxury,  and  the  values  stored  in  most  hold  circuits  must  be  retained  as  state  variables. 

Consider  a  sample/hold  circuit  with  period  T\  =  pit,  input  u,  and  output  y.  Let  x(nt)  be 
the  value  stored  in  the  hold  circuit  at  time  t  =  nt.  For  times  t  =  not  at  which  sampling  occurs,  the 
next  state  is  updated;  as  a  consequence  of  Assumption  2.1(d),  the  output  is  set  to  the  input: 


(a,b,c,d)(nt)  = 


x((no+l)x)  =  u(not)  and  y(nox)  =  u(nox). 

At  all  other  times,  the  state  remains  unchanged  and  the  output  is  set  to  the  value  of  the  state: 

x((n+l)x)  =  x(nt)  and  y(nx)  =  x(nx). 

Thus,  for  a  sample/hold  with  period  Ti, 

x((n+l)x)  =  a(nx)  x(nx)  +  b(nx)  u(nx) 
y(nx)  =  c(nx)  x(nx)  +  d(nx)  u(nx), 

where 

(0,1, 0,1)  ,  if  n  mod  pi  =  0 
(1,0, 1,0)  ,  if  n  mod  pi  *  0 

Note  that  if  pi  =  1,  the  representation  becomes  time  invariant,  with  corresponding  state  and  output 
equations  x((n+l)x)  =  u(nx),  y(nx)  =  u(nx).  Clearly,  x(nx)  is  superfluous  in  this  case,  and  states 
need  not  be  assigned  to  sample/holds  of  period  x. 

The  following  convention  is  explicidy  stated  for  clarity. 

Convention  3.1:  (Intermediate  state  values  of  discrete-time  systems) 

Given  a  single-rate  discrete-time  subsystem  with  state  x  and  state  transitions 
specified  at  times  kTj  or  a  multirate  discrete-time  subsystem  with  state  x  specified 
in  Trexpanded  form,  define 

x(t)  =  x(kTi),  kTl^t<((k+l)T1). 

♦  ♦  ♦ 

Due  to  the  multirate  nature  of  the  surrounding  elements,  a  T! -varying  discrete-time 
subsystem  must  sample  its  inputs  at  times  t  =  kTj  and  maintain  its  output  at  y(t)  =  y(kTj)  for 
kTj  <  t  <  (k+l)Ti.  Since  u(kTi)  is  required  at  time  t  =  kTi+(pj-l)x  to  determine  the  state 
transition,  a  state  variable  must  be  defined  to  retain  u(kTj).  Let  a  periodic  Ti-varying  single-rate 
discrete-time  subsystem  with  realization  (A(kTi),B(kTi),C(kTi),D(kTi))  have  state  x  €  Rn,  input 
uelm  and  output  y(kTi).  Define  the  augmented  state 

x((kpi+i)x) 


z((kpt+i)x)  = 


,  0  <  i  <  pi,  k  =  0,1,... 


u(kT0  J 

A  x-varying  representation  of  this  system  which  satisfies  the  state  and  output  equations 

z(t+x)  =  a(t)  z(t)  +  b(t)  u(t) 
y(t)  =  c(t)  z(t)  +  d(t)  u(t) 

for  times  t  =  (kpi-H)x,  0  <  i  £  pi  and  k  =  0,1,...,  is 
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(a,b,c,d)((kpi+i)t)  =  ■ 


(  T  *n  0  1  r 0  I  > 

,  [COcTO  0],  D(kTO  L  if  i  =  0 
VL”  UJ  L^mJ  ) 

(  r  i„  0  i  roi  A 

_  _  »  ,  [C(kTj)  D(kT!)],  0  ,  if  0  <  i  <  Pl - 1 

v  L  o  im  J  LOJ 

(  T  A(kTO  B(kTi)  "I  rOl 

,  [C(kT,)  D(kT,)],  0  Lifi  =p, 

V  L  0  ) 


=  Pi  -  1 


(3.1) 

The  reader  may  wish  to  trace  through  the  above  representation  and  note  the  following  facts. 
u(kTi)  is  updated  at  the  time  it  is  available,  i  =  0,  and  the  following  output  values  and  state 
transition  are  determined  from  this  stored  value.  The  output  is  held  at  y((kpi+i)x)  =  y(kTj)  for 
0  <  i  <  pi .  This  x-varying  representation  is  equivalent  to  the  Ti-varying  representation  in  the 
sense  that  for  any  z(0)  e  Rn+m  and  input  sequence  { u(jx):  j  =  0,1,...},  the  resulting  output 
sequence  {y(jx)}  and  the  x  portion  of  the  state  trajectory  { z(jx) }  of  the  x-varying  representation  are 
identical  to  samples  at  times  t  =  jx  of  the  output  sequence  { y(kT i ) }  and  the  state  trajectory 
{x(kTi)},  where  x(jx)  is  obtained  from  x(kTi)  via  Convention  3.1,  of  the  Tj-varying 
representation  with  initial  state  x(0)  and  input  u  =  { u(jx) } . 

A  periodic  multirate  discrete-time  system  specified  in  Ti -expanded  form. 


x(kTi)  *  A(kTi)  x(kTi)  +  B(kT0  U(kTj) 

Y(kTj)  =  C(kTi)  x(kT0  +  D(kT0  U(kT0, 

may  also  be  expressed  as  a  x-varying  system.  It  is  necessary  that  D(kTi)  be  structured  so  that  this 
system  is  causal.  The  state  of  the  system  must  be  augmented  by  U(kTi).  Although  the  specific 
details  of  a  x-varying  representation  for  a  general  multirate  discrete-time  system  specified  in  time 
expanded  form  are  too  involved  to  present  here,  the  basic  philosophy  is  the  same  as  that  for  a 
single-rate  system.  The  values  comprising  U(kTi)  are  loaded  into  the  augmented  state  as  they 
become  available,  and  the  state  transition  occurs  during  the  last  x  interval  of  each  Tj  interval,  as 
specified  in  Convention  3.1.  Strictly  as  a  matter  of  convenience  for  developments  in  Section  3.3. 
the  following  convention  should  be  satisfied  by  a  x-varying  representation  of  a  multirate  discrete¬ 
time  system  specified  in  time  expanded  form. 

Convention  3.2: 


Portions  of  the  augmented  state  corresponding  to  input  values  that  are  not 
yet  available  are  assigned  the  value  zero. 


♦  ♦  ♦ 


Example  3.1: 

A  SISO  multirate  system  with  a  single  state  is  specified  in  6x-expanded 


form: 


fi5 


a 


x((k+l)6t)  =  a  x(k6t)  +  [bj  l*]  U(k6t) 

U(k6x). 


"Cl  ' 

"  d, 

o  - 

Y(k6x)  = 

c2 

x(k6x) + 

d2 

0 

_C3  _ 

_  d3 

da  _ 

The  input  u  is  sampled  with  period  3x  and  the  output  y  has  period  2x: 

U(k6x)  = 


_  _ 

y(k6x) 

u(k6x) 

,  Y(k6x)  = 

y(k6x+2x) 

u(k6x+3x) 

_y(k6x+4x)_ 

Note  the  two  0's  in  the  D  matrix  required  for  causality.  Define  the  3  x  1  augmented 
state  z: 

rx((k6+i)x) 


z((k6+i)x)  = 


,  0  <  i  <  6,  k  =  0,1 . 


U(k6x) 

As  this  system  has  period  6x,  six  sets  of  matrices  are  required  for  the  x-varying 
representation.  These  matrices  are  given  as  a(t),  b(t),  c(t),  and  d(t)  for  each  t 
below. 


-10  0- 

-0- 

t  =  k6x : 

0  0  0 

r 

1 

.0  0  0. 

_0_ 

,  [cj  0  0],  dj 


-10  0- 

-o- 

t  =  k6x  +  x  : 

0  1  0 

♦ 

0 

.0  0  0. 

.0. 

-10  0- 

-0- 

t  =  k6x  +  2x : 

0  1  0 

* 

0 

.0  0  0  . 

.0. 

-10  0- 

-0- 

t  =  k6x  +  3x  : 

0  1  0 

» 

0 

.0  0  0  . 

_  1  _ 

-10  0- 

-o- 

t  =  k6x  +  4x  : 

0  1  0 

* 

0 

.001. 

.0. 

,  [cj  d!  0],  0 


,  [c2d20],  0 


,  [c2  d2  0],  0 


,  [C3  d3  da],  0 
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A  continuous-time  subsystem  whose  inputs  and  outputs  are  sampled  and  held  at  multiple 
rates  is  readily  converted  to  x-varying  form.  The  continuous-time  subsystem  is  merely  discretized 
with  period  x  and  the  surrounding  sample/hold  devices  govern  the  flow  of  input  and  output  values 
to  and  from  the  multirate  sampled  subsystem.  The  following  procedure  details  the  calculation  of 
the  x-varying  representation  of  a  multirate  discrete-time  system. 

Procedure  3.1: 

a.  Verify  that  Assumption  2.1  is  met,  with  the  exception  that  periodic 
single-rate  and  multirate  discrete-time  subsystems  may  be  present.  Partition  the 
entire  system  into  subsystems  which  satisfy  Assumption  2.1(e). 

b.  Use  Procedure  2.1  to  normalize  all  sampling  periods  and  periods  of 
periodic  subsystems  and  find  T  =  Px. 

c.  Find  a  state  space  realization  for  each  continuous-time  subsystem  and 
single-rate  discrete-time  subsystem. 

d.  Discretize  with  period  x  each  continuous-time  subsystem  and  find  a 
x-varying  representation  for  each  sample/hold  device  and  discrete-time  subsystem  in 
accordance  with  Convention  3.1  and  Convention  3.2  as  described  in  this  section. 

e.  For  each  j  g  {0,1,...,P-1 },  associate  the  t  =  jx  value  of  the  x-varying 
representation  of  each  subsystem  and  sample/hold  device  with  the  corresponding 
subsystems  and  sample/hold  devices  in  the  block  diagram  of  the  system.  Reduce 
this  block  diagram  to  a  single  state  space  equation  as  if  it  were  the  block  diagram  of 
a  single-rate  system.  This  gives  a  T-periodic,  x-varying  representation  for  the  entire 
multirate  system:  (A(nx),B(nx),C(nx),D(nx)). 

♦  ♦  ♦ 

As  with  Procedure  2.4,  difficulties  may  arise  when  performing  Procedure  3.1(e). 

In  certain  instances,  Procedure  3.1(d)  can  be  modified  in  order  to  reduce  the  number  of 
states  of  the  x-varying  representation.  For  example,  if  a  signal  is  connected  to  a  sample/hold  with 
period  Tj  or  to  a  single-rate  discrete-time  subsystem  which  is  T\ -sampling  the  signal  and  this 
signal  is  the  output  of  a  sample/hold  with  period  T2,  a  Tj-held  output  of  a  discrete-time  system,  or 
the  sum  of  such  signals  and  T2/T1  g  N,  then  a  state  variable  need  not  be  assigned  to  the 

sample/hold  with  period  Tj  or  to  retain  the  Tj-sampled  input  value  of  the  discrete-time  subsystem. 
In  this  case,  a  x-varying  representation  for  a  Ti-varying  single-rate  discrete-time  subsystem  is,  in 


3? 


the  notation  of  equation  (3.1), 


i 

I 

f 


(a,b,c,d)((kpi+i)x) 


(In,  0,  C(kTi),  D(kTj))  ,  if  0  <  i  <  Pl- 1 
(A(kTj),  B(kTj),  C(kTi),  D(kTO)  ,  if  i  =  Pl- 1 


(3.2) 


If  some  of  the  inputs  to  the  discrete-time  subsystem  have  this  property  and  others  do  not,  the 
proper  t-varying  representation  is  a  hybrid  of  (3.1)  and  (3.2).  The  justification  for  this 
modification  is  that  when  T2/T1  e  N,  the  inputs  to  the  sample/hold  of  period  Ti  or  the  Ti -sampled 
inputs  of  a  discrete-time  subsystem  are  held  constant  over  the  times 


t  =  kTi,  kTi  +  x,...,  kTi  +  (pi-l)x; 


u(kTi)  is  available,  if  necessary,  at  times  after  t  =  kTi  for  determining  outputs  and  state  transitions. 

A  slight  modification  of  Procedure  3.1  produces  a  t/q-varying  representation,  qe  N.  By 
introducing  a  fictitious  sampling  period  To  =  x/q  into  the  normalization  process  in  Procedure 
3.1(b),  a  different  set  of  parameters  results:  x'  =  t/q,  P'  =  qP,  and  p/  =  qp;.  Utilizing  these 
parameters  in  place  of  x,  P,  and  p,  in  steps  (c),  (d),  and  (e)  of  Procedure  3.1  then  gives  a 
x'-varving  representation. 


3.2  The  mx-varying  Representation 

Obtaining  an  mx-vaiying  representation,  m  €  N,  for  a  multirate  system  by  direct  means  can 
be  somewhat  tedious.  However,  a  method  will  be  outlined  whereby  an  mx-varying  representation 
can  be  obtained  via  the  x-varying  representation  in  a  straightforward  manner.  In  the  following 
procedure,  assume  that  the  system  is  single-bundle  input,  single-bundle  output. 

Procedure  3.2: 

a.  Apply  Procedure  3.1  to  obtain  (A(nx),B(nx),C(nx),D(nx)). 

b.  Let  r  =  LCM  {m,P}/m.  For  k  =  0,1,...,  calculate 

a(kmx)  =  A((km+m-l)x)...A((km+l)x)A(kmx), 

P(kmx)  =  [A((km+m-l)x)...A((krn+l)x)B(kmx)  I ... 

...  I  A((km+m-l)x)B((km+m-2)x)  I  B((km+m-l)x)], 


Y(kmx)  = 


C(kmx) 

C((km+l)x)A(kmx) 


|_C((km+m- 1  )x)A((km+m-2)x)...  A(kmx)  J 


and  the  block  m  x  m  matrix  8(kmx),  where 


Sjj(kmT)  - 


55  $ 


Vl» 

’V 

•M 

w 


lV  •' 

s  ' 


t 


ft  ft 


r  0  ,  if  i  <  j 
D((km+i-l)x)  ,  if  i  =  j 

8ij(kmx)  =  < 

C((km+i-l)x)B((km+j-l)x)  ,  if  i  =  j  +  1 
t  C((km+i-l)x)A((km+i-2)x)...A((km+j)x)B((km+j- 1  )x)  ,  if  i  >  j  +  1 

Since  a(.),  (3(.),  y(.),  and  8(.)  are  periodic  with  period  rmx,  only  values  for 
k  =  0,...,  r  -  1  are  required.  With  z(nx),  u(nx),  and  y(nx)  as  the  state,  input,  and 
output  of  the  x-varying  representation,  the  mx-varying  representation  satisfies  the 
following  state  and  output  equation  for  k  =  0,1,... : 

z((k+l)mx)  =  a(kmx)  z(kmx)  +  (3(kmx)  U(kmx) 

Y(kmx)  =  y(kmx)  z(kmx)  +  5(kmx)  U(kmx). 

Y(kmx)  and  U(kmx)  are  the  mx-expanded  versions  of  u(nx)  and  y(nx). 

♦  ♦  ♦ 

Iteration  of  the  equations  governing  the  x-varying  representation  confirms  the  validity  of 
Procedure  3.2.  Note  the  similarity  between  these  formulas  and  those  for  the  time  expansion  of  a 
periodic  discrete-time  system  in  Section  2.7.  Y(kmx)  and  U(kmx)  in  Procedure  3.2  should  not  be 
confused  with  the  mx-expanded  versions  of  the  output  and  input  of  the  original  system.  Y(kmi) 
and  U(kmx)  are  mx-expanded  versions  of  x-sampled  versions  of  the  original  system's  input  and 
output.  In  fact,  the  representation  produced  by  Procedure  3.2  is  the  mx-varying  representation  of 
the  multirate  system  with  samplers  of  period  X  inserted  at  each  input  and  output  of  the  original 
system.  These  added  samplers  do  not  alter  the  behavior  of  the  system.  Since  the  inputs  and 
outputs  of  a  multirate  system  are  generally  sampled  with  periods  greater  than  x,  many  of  the 
components  of  U(kmx)  are  not  accessed  by  the  mx-varying  representation  and  many  of  the 
components  of  Y(kmx)  are  duplicates.  This  redundancy  serves  to  pad  the  inputs  and  outputs  of  the 
mx-varying  representation  to  a  size  which  is  uniform  in  time. 

Example  3.2: 

Consider  finding  the  3x-varying  representation  of  the  multirate  system  in 
Example  3.1.  Since  r  =  6/3  =  2.  only  two  values  for  each  of  a.  (3,  y,  and  5  are 
needed.  The  x-varying  representation  of  this  system  was  determined  in  Example 
3.1;  Procedure  3.2(a)  is  finished.  By  straightforward  calculations, 
k  =  0: 

■10  0*1  ro  o  oi  rci00-|  rd,  o  fl- 

a(0)  =  0  0  0  ,  (3(0)  =  100  ,7(0)=  c,  0  0  .5(0)=  d,  0  0  . 

_000j  L  0  0  0  J  L  c2  0  0  J  [_  d2  0  0  _ 
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Y(kmt)  =  y(t)  ,  and  U(kmx)  =  u(x) 
_y(2x)J  l“(2x). 


k=  1: 

“  a  bi  0  “I  r  b2  0  0  “I  r*c2  d2  0  "I  r  0  0  O' 

a(3x)  =010,  p(3t)  =  0  0  0,  y(3x)  =  c3  d3  0  ,  5(3x)  =  d4  0  0 

_0  0  Oj  L  1  0  0 J  [_  c3  d3  0  J  1_ d4  0  0 _ 


Thus, 


y(3x)  ~u(3x) 

Y(kmx)  =  y(4x)  ,  and  U(kmx)  =  u(4x) 
„y(5x)_  _u(5x) 


(a,p,v,8)  (k(3x))  = 


(a.p.y.S)  (0)  ,  if  k  =  0,2,4,... 
(a,p,Y,5)  (3x)  ,  if  k  =  1,3,5,... 


Since  y  has  period  2x,  it  should  be  the  case  that  y(0)  =  y(x),  y(2x)  =  y(3x), 
and  y(4x)  =  y(5x),  which  can  be  verified.  Note  that  only  u(0)  and  u(3x)  are 
utilized,  which  is  consistent  with  the  fact  that  u  is  updated  with  period  3x  by  the 
original  multirate  system. 


Procedure  3.2  determines  the  mx-varying  representation  of  a  single-bundle  input,  single¬ 
bundle  output,  multirate  system.  The  mx-varying  representation  of  a  multirate  system  with  inputs 
and  outputs  which  are  vectors  of  bundles  results  from  partitioning  B(nx),  C(nx),  and  D(nx)  and 
repeated  application  of  Procedure  3.2(b)  in  a  manner  analogous  to  the  time  expansion  of  such 
multirate  systems.  As  noted  earlier,  x/q-varying  representations,  qe  N,  may  be  obtained  using 
Procedure  3.1.  It  follows  from  the  developments  in  this  section  that  if  r  e  Q  and  r  >  0,  Procedure 
3.2  may  be  used  to  calculate  the  rx-varying  representation  of  a  multirate  system,  since  r  =  m/q  for 
some  integers  m  and  q. 


3.3  Corresponding  Representations 

A  trait  common  to  all  periodic  representations  of  multirate  systems,  but  not  the  T-expanded 
representation,  is  that  states  corresponding  to  the  contents  of  some  of  the  hold  circuits  in  the 
system  must  be  included  in  the  composite  state  vector  of  the  system.  These  added  states  are  of  a 


nondynamic  nature,  and  the  reader  may  correctly  predict  that  the  lower-dimensional  T-expanded 
representation  of  a  given  multirate  system  shares  many  of  the  qualitative  properties  of  a  periodic 
representation  of  that  system.  The  representations  described  in  this  section  serve  as  a  crucial  step 
in  the  investigation  in  Chapter  4  of  the  manner  in  which  the  T-expanded  representation  relates  to 
periodic  representations.  The  properties  and  results  discussed  in  this  section  and  Chapter  4 
regarding  T-expanded  and  x-varying  representations  could  be  adapted  to  almost  any  penodic 
representation  of  the  types  of  multirate  systems  under  consideration. 

For  a  given  multirate  system,  a  relationship  exists  between  the  matrices  A(nx)  of  the 
x-varying  representation  obtained  from  Procedure  3.1  and  the  matrix  Ae  of  the  T-expanded 
representation  obtained  from  Procedure  2.4,  provided  that  two  particular  steps  of  these  procedures 
are  performed  similarly.  To  ensure  that  such  representations  exist,  let  the  following  assumption  be 
satisfied. 

Assumption  3.1:  (well-posedness) 

The  calculations  necessary  to  complete  Procedure  2.4(i)  and  Procedure 
3.1(e)  can  be  performed  for  the  multirate  systems  under  consideration  and  involve 
matrices  with  bounded  elements. 

♦  ♦  ♦ 

The  state  of  the  T-expanded  representation,  x(kT),  is  a  composite  of  the  states  of  all  the  dynamic 
subsystems  in  the  multirate  system.  The  state  of  the  x-varying  representation,  z(nx),  is  a  composite 
of  the  states  of  all  the  dynamic  subsystems  in  the  system  and  hold  states  (states  assigned  to 
sample/hold  circuits  and  states  introduced  to  retain  input  values  of  discrete-time  subsystems.)  A 
T-expanded  and  a  x-varying  representation  of  a  multirate  system  will  be  called  corresponding 
representations  if  these  representations  are  obtained  from  Procedure  2.4  and  Procedure  3.1, 
respectively,  and  the  following  convention  is  satisfied. 

Convention  3.3: 

a.  The  same  system  partition  and  the  same  state  space  realization  for  each 
subsystem  is  used  in  step  (c)  of  Procedures  2.4  and  3.1. 

b.  When  reducing  the  composite  system  to  a  single  state  space  equation  in 
Procedure  3.1(e),  z(nx)  is  partitioned  into  dynamic  states,  w(nx),  and  hold  states, 
h(nx): 

w(nx) 

z(nx)  = 

.h(nx)_ 

In  addition.  Procedure  3.1(e)  and  Procedure  2.4(i)  are  performed  so  that  the  state  of 
any  given  subsystem  occupies  the  same  components  of  x  and  w. 

♦  ♦  ♦ 
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If  corresponding  representations  have  initial  conditions  z(0)  and  x(0),  where  x(0)  =  w(0), 
and  are  subjected  to  the  same  input,  then  x(kT)  =  w(kPx),  k  =  0,1,...,  because  x  and  w  each 
represent  the  dynamic  states  of  the  system  and  are  structured  identically.  Hereafter,  this  fact  is 
emphasized  by  writing 


z(nx)  = 


x(nt) 


Lh(nx)J 

where  x(nx)  is  interpreted  through  Convention  3.1  as  representing  the  dynamic  states  of  the 
subsystems  in  the  multirate  system  at  time  t  =  nx. 

Let  x  be  5  x  1  and  h  be  rj  x  1.  As  a  consequence  of  the  manner  in  which  x-varying 
representations  were  defined  in  Section  3.1,  at  time  t  =  0  each  subsystem  and  sample/hold 
determines  its  next  state  and  current  output  solely  from  the  current  value  of  its  dynamic  states  (if 
any)  and  inputs.  Thus,  at  t  =  0  the  next  state  z(x)  is  independent  of  every  component  of  h(0).  If 
the  inputs  to  the  system  are  zero,  from  z(x)  =  A(0)z(0),  it  follows  that  for  some  Ei  e  R5x5  and 
E2  €  RnxS, 


A(0)  = 


‘E!  0‘ 
_E2  0. 


Expressing  z(Px)  =  z(T)  for  arbitrary  z(0)  as 


z(px)  =  A((P-l)x)...A(x)A(0)z(0) 


'G  O' 
H  0 


z(0) 


for  some  G  e  R5x5  and  H  e  R^x5,  it  follows  that  x(Px)  =  Gx(0)  for  arbitrary  x(0).  For  the 
corresponding  T-expanded  representation,  x(T)  =  Aex(0)  and  x(kT)  =  x(kPx).  Thus, 


Aex(0)  =  x(T)  =  x(Px)  =  Gx(0) 


for  arbitrary  x(0),  which  implies  that  G  =  Ae.  The  following  theorem  results  from  the  discussion 
above. 

Theorem  3.1: 

For  corresponding  T-expanded  and  x-varying  representations  of  a  multirate 
system  and  some  Ej  e  R5xS  and  E2,  H  6  R1^5, 


A(0)  = 


'E!  O' 
_E2  0. 


and  A((P-l)x)...A(x)A(0) 


’  Ae  O' 
H  0 


♦  ♦  ♦ 


A  T-expanded  and  a  x-varying  representation  of  a  multirate  system  will  be  called 
completely  corresponding  representations  if  these  representations  are  corresponding 
representations  and  the  following  convention  is  satisfied. 
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Convention  3.4: 

a.  Samplers  of  period  x  are  placed  on  each  input  and  output  of  the  block 
diagram  of  the  system  before  finding  the  T-expanded  representation. 

b.  Once  the  T-expanded  representation  has  been  obtained,  its  input  and 
output,  U(kT)  and  Y(kT),  are  restructured  to  form  a  new  input  and  output,  U(kPx) 
and  Y(kPx),  given  by 


'  u(kPx)  ' 

y(kPx) 

U(kPx)  = 

u(kPx+x) 

and  Y(kPx)  = 

y(kPx+x) 

_u((k+l)Px-x)_ 

_y((k+l)Px-x)_ 

In  essence,  the  input  and  output  of  the  system  are  each  treated  as  a  single  bundle  of 
signals.  During  this  operation,  the  rows,  columns,  and  elements  of  the  matrices  of 
the  representation  must  be  permuted  according  to  the  reordering  of  the  inputs  and 
outputs. 

♦  ♦  ♦ 

Let  the  completely  corresponding  T-expanded  representation  be  denoted  by 

x((k+l)T)  =  Ae  x(kT)  +  Bcc  U(kPx) 

Y(kPx)  =  Ccc  x(kT)  +  Dcc  U(kPx) 

The  primary  difference  between  a  corresponding  and  a  completely  corresponding  T-expanded 
representation  is  Convention  3.4(a);  however,  the  added  samplers  do  not  alter  the  behavior  of  the 
system.  In  fact,  the  columns  of  Bcc  are  permutations  of  the  columns  of  Be  interspersed  with 
columns  of  zeroes,  and  the  rows  of  Qc  contain  duplicates  and  are  drawn  from  the  rows  of  Ce.  For 
a  given  multirate  system,  a  t-varying  representation  which  is  a  corresponding  representation  and 
one  which  is  a  completely  corresponding  representation  are  identical. 

The  manner  in  which  completely  corresponding  representations  are  related  is  revealed  by 
defining  the  following  quantities  from  the  matrices  (A(nx),B(nx),C(nx),D(nt))  of  the  T-varying 
representation: 

a  =  A((P-l)x)...A(x)A(0), 

P  =  [A((P-l)x)...A(x)B(0) ! ...  I  A((P-1)t)B((P-2)t)  I  B((P-l)t)l, 


C(0) 

C(T)A(0) 


C((P-1)t)A((P-2)t)...A(0) 


and  the  block  P  x  P  matrix  8,  where 
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5ijH 

J  C((i-l)i)B((j-l)x)  ,  if  i  =  j  +  1 

.  C((i-l)t)A((i-2)t)...A(jx)B((j-l)t)  .  if  i  >  j  +  1 
Note  the  similarity  of  these  quantities  to  the  formulas  in  Procedure  3.2  for  a  Px-varying 
representation  of  a  multirate  system  whose  inputs  and  outputs  are  treated  as  a  single  bundle.  It 
follows  that 

z((k+l)Px)  =  a  z(kPx)  +  (3  U(kPt) 

Y(kPx)  =  y  z(kPx)  +  5  U(kPx). 

For  completely  corresponding  representations  with  dynamic  states  which  are  identical  at  t  =  0  and 
input  U(kPt),  not  only  will  x(kT)  =  x(kPx),  as  with  corresponding  representations,  but  in  addition 
their  outputs  Y(kPx)  will  be  identical.  In  particular,  with  z(0)  =  0,  x(0)  =  0,  and  U(0)  arbitrary. 


x(T)  =  Bcc  U(0),  Y(0)  =  Dcc  U(0) 


and 


z(Px)  = 


x(Pt) 

h(Px) 


=  p  U(0),  Y(0)  =  5  U(0). 


Since  U(0)  is  arbitrary  and  x(T)  =  x(Pt),  P  must  be  of  the  form 


for  some  matrix  J  and  5  =  Dcc.  By  the  manner  in  which  x-varving  representations  were  defined  in 
Section  3.1,  at  time  t  =  0  each  subsystem  and  sample/hold  determines  its  next  state  and  current 
output  solely  from  the  current  value  of  its  dynamic  states  (if  any)  and  inputs.  Thus,  for  some 
matrices  Ej,  Ei,  and  E3, 

"  E  0 " 

A(0)  =  1  and  C(0)  =  [E3  0]. 

_  E2  0 

By  definition  of  y,  it  follows  that  for  some  matrix  E4,  y  =  [E4  0J.  Consider  setting  U(0)  =  0  and 
letting  z(0)  and  x(0)  be  arbitral.  Then 

rx(0)i 

yz(0)-[E40l[  J  =  Y(0)-Cccx(0v. 

therefore,  y  =  rCcc  0].  The  following  theorem  summarizes  the  a  .  a  c  di^  ussion. 

Theorem  3.2: 

For  completely  corresponding  T-expanded  and  t  \  .iryinc  representations  of 
a  multirate  system  and  some  matrices  H  and  J, 


CHAPTER  4 


ANALYSIS  OF  MULTIRATE  SYSTEMS 

The  representations  developed  in  the  previous  two  chapters  provide  a  convenient  means  of 
assessing  important  qualitative  characteristics  of  multirate  systems.  Of  particular  value  is  the  fact 
that  the  stability,  controllability,  reconstructibility,  stabilizabilitv,  and  detectability  of  the  i-varying 
representation  may  be  determined  by  examining  the  corresponding  T-expanded  representation  of 
the  system. 

In  this  chapter,  only  systems  that  satisfy  Assumption  3.1  (well-posedness)  are  considered. 
Throughout,  denote  the  x-varying  representation  of  the  multirate  system  by 

z((n+l)x)  =  A(nx)  z(nx)  +  B(nt)  u(nx) 
y(nx)  =  C(nx)  z(nx)  +  D(nx)  u(nx), 

where  (A,B,C,D)  (nx)  is  T-periodic,  and 

x(nx) 

z(nx)  =  €  . 

.h(nx)_ 

Denote  the  corresponding  T-expanded  representation  of  this  system  (see  Section  3.3  for  details  i  by 

x((k+l)T)  =  Ae  x(kT)  +  Be  U(kT) 

Y(kT)  =  Ce  x(kT)  +  De  U(kT) 

and  the  completely  corresponding  T-expanded  representation  of  this  system  by 

x((k+l)T)  =  Ae  x(kT)  +  Bcc  L'(kPx) 

Y(kPx)  =  Ccc  x(kT)  +  Dcc  U(kPx). 

where  x  €  ?.5  and  T  =  Px.  By  Assumption  3.1,  the  elements  of  (A.B.C.D)  (nx),  ( Ae.Be.Cc,De). 
and  (Ae,Bcc,Ccc,Dcc)  are  bounded  and  well-defined. 

4.1  Stability 

Clearly,  a  necessary  condition  for  uniform  asymptotic  stability  of  a  multirate  system  is  that 
its  T-expanded  representation  be  asymptotically  stable,  i.e..  all  the  eigenvalues  of  A0  have 


magnitudes  strictly  less  than  one.  It  turns  out  that  this  is  also  a  sufficient  condition  for  uniform 
asymptotic  stability  of  the  corresponding  x-varying  representation. 

Theorem  4.1: 

If  the  eigenvalues  of  Ae  have  magnitudes  less  than  1,  then  the  x-varying 
representation  is  uniformly  asymptotically  stable. 

Proof: 

Let  Z(n,no,zo)  be  the  zero  input  state  trajectory  at  time  t  =  nx  of  the 
x-varying  representation  resulting  from  the  initial  state  zo  at  time  t  =  nox.  It  suffices 
to  show  that  V  no  e  N  and  zo  e  Xs+Tl,  Hm  Z(n,no>zo)  =  0,  as  this  result  and  the 

n — >oo 

periodicity  of  the  x-varying  representation  establish  uniformity. 

Assume  the  eigenvalues  of  Ae  have  magnitudes  less  than  1.  Define  the 
monodromy  matrix  at  time  nx  for  the  x-varying  representation  as 

O(n)  =  A((n+P-l)x)...A((n+l)x)A(nx). 

Note  that  Otn+kP)  =  O(n)  for  k  =  0,1,...  and  that  under  zero  input 
z((n+P)x)  =  0(n)z(nx)  for  n  =  0,1,...  . 

By  Theorem  3.1,  for  some  H  6  R’l*5, 

Ae  0 

H  0 

Thus,  each  eigenvalue  of  0(0)  has  magnitude  less  than  one.  Recall  that  for  square 
matrices  V  and  W,  VW  and  WV  have  the  same  eigenvalues.  Repeated  application 
of  this  fact  leads  to  the  conclusion  that  0(0),  0(1),...,  and  O(P-l)  all  have  identical 
eigenvalues.  Let  i  e  {1 . P)  and  k  =  0,1,...  and  note  that 

Z(kP-t-i+no,n0,zo)  =  [O(no+i)]k  A((n0+i-i)x)...A(n0x)zo. 

For  any  zo  e  a5+ti  and  any  no  e  N,  A((no4-i-i)x)...A(noX)zo  is  bounded  for  each 

i  e  { 1 . P}.  Since  the  eigenvalues  of  each  O(n)  have  magnitudes  less  than  1.  for 

each  i  €  { 1,...,P)  Z(kP+i+no,no,zo)  is  bounded  and  iim  Z(kP+i+no,no,zo)  =  0; 

K — *  oo 


0(0)  = 


Zfn.no,zo)  is  partitioned  into  P  bounded  convergent  subsequences.  Thus, 

Z(n,no,zo)  is  bounded  and  for  any  i  e  { 1 . P)  and  £  >  0,  there  is  K(i.e)  e  N 

such  that  ||  Z(kP+i+no,no.zo)  II  <  £  V  k  >  K(i,£).  It  follows  that  for  any  £  >  0  there 

is  K(e)  =  1  +  max{K(i,£):  i  e  { 1 P} }  such  that  II  Z(n.no.zo)  II  <  £  for  all 

n  >  K(£)P  +  no,  which  implies  that  lim  Z(n,no,zo)  =  0. 


The  utility  of  Theorem  4.1  stems  from  the  fact  that  T-expanded  representations  are 
generally  easier  to  calculate  and  are  of  lower  order  than  t-varying  representations.  Theorem  4. 1 
also  insures  that  x(nt),  rather  than  merely  x(kT),  decays  asymptotically  to  zero  for  any  initial 
condition  if  the  T-expanded  representation  is  asymptotically  stable.  In  fact,  consider  collecting  the 
states  of  all  continuous-time  subsystems  in  the  multirate  system  into  the  composite  state  xc(t).  For 
k  =  0,1,...  and  0  <  to  <  Pt,  write 

xc  (kPx+to)  =  ^(to)  z(kPx),  (4.1) 

where  T^to)  is  a  dim(Xc)  x  (S+rj)  matrix.  Under  the  assumption  that  llvF(to)ll  is  bounded  for  each 
to  s  [0,Pt),  if  the  eigenvalues  of  Ae  have  magnitude  less  than  1,  it  follows  from  (4.1)  and 
Theorem  4. 1  that  as  t  — »  «>,  xc(t)  — » 0.  As  a  consequence  of  Theorem  4. 1  and  the  satisfaction  of 
Assumption  3.1  by  the  systems  under  consideration,  if  the  eigenvalues  of  A«  have  magnitudes  less 
than  1,  the  t-varying  and  T-expanded  representations  are  BIBO  stable  (Chen,  1987). 

4.2  Controllability  and  Reconstructibility 

The  controllability,  reachability,  reconstructibility,  and  observability  of  linear  periodic 
discrete-time  systems  have  been  investigated  in  several  recent  works;  see  Grasselli  (1984),  Bittani 
and  Bolzern  (1985),  Bittani  and  Colaneri  (1986),  Bittani  and  Guardabassi  (1986)  and  the 
references  cited  therein.  The  t-varying  representation  provides  a  T-periodic  description  of  a 
multirate  system  amenable  to  analysis  by  the  methods  developed  in  these  references.  Such  an 
analysis  may  be  quite  involved  because  the  properties  of  controllability,  reachability,  etc.  of  a  linear 
periodic  system  are  in  general  enjoyed  by  only  a  time-varying  subspace  of  the  entire  state  space. 
The  analysis  here  presents  two  instances  in  which  the  t-varying  representation  inherits  a  property 
of  its  corresponding  T-expanded  representation. 

The  informal  definitions  which  follow  are  consistent  with  those  of  Grasselli  (1984).  A 
linear,  discrete-time  system  is  said  to  be  controllable  if  there  exist  inputs  which  drive  the  system 
from  any  given  initial  state  at  any  given  time  to  the  zero  state  in  a  finite  amount  of  time.  A  linear, 
discrete-time  system  is  said  to  be  reachable  if  there  exist  inputs  which  drive  the  system  to  any 
given  terminal  state  at  any  given  point  in  time  from  the  zero  state  in  a  finite  amount  of  time.  A 
linear,  To-varying  discrete-time  system  is  said  to  be  reconstructive  if  the  state  of  the  system  at 
any  time  t  =  mTo  can  be  determined  from  knowledge  of  the  inputs  and  outputs  of  the  system  over 
some  finite  interval  of  time  prior  to  t  =  mTo.  A  linear.  To-varying  discrete-time  system  is  said  to  be 
observable  if  the  state  of  the  system  at  any  time  t  =  mTo  can  be  determined  from  knowledge  of 
the  inputs  and  outputs  of  the  system  over  some  finite  interval  of  time  beginning  at  t  =  mTy. 


Determination  of  the  controllability,  reachability,  reconstructibility,  or  observability  of  the 
T-expanded  representation  can  be  accomplished  by  considering  (Ae,Be,Ce)  as  a  time-invariant 
triple.  A  consequence  of  the  Cayley-Hamilton  theorem  is  that  at  most  5  T-expanded  inputs  are 
required  to  drive  the  state  x(kT)  to  zero  if  the  T-expanded  representation  is  controllable  or  to  the 
desired  terminal  state  if  the  T-expanded  representation  is  reachable.  In  addition,  at  most  8 
T-expanded  outputs  are  required  to  determine  a  present  state  x(kT)  of  the  T-expanded 
representation  if  it  is  reconstructible  or  to  determine  an  initial  state  x(kT)  of  the  T-expanded 
representation  if  it  is  observable. 

Theorem  4.2: 

The  x-varying  representation  is  controllable  if  the  corresponding  T-expanded 
representation  is  controllable,  and  any  initial  condition  z(nox)  =  zo  e  RTi+5  can  be 
driven  to  the  zero  state  by  time  t  =  not  +  (S+1)T, 

Proof: 

Let  no  =  koP  -  i,  i  e  {0,1,...,P-1 }.  If  i  >  0  and  the  inputs  at  times 
not,...,  (no+i-l)t  are  zero,  the  system  progresses  from  zo  to  some  state 

r  x(k0pt)i 


z(koPt)  = 

_h(k0Pt)_ 

where  x(koPt)  =  x(koT)  e  Rs.  Since  the  T-expanded  representation  is  con¬ 
trollable,  there  is  an  input  sequence  which  drives  x(koT)  to  the  origin  by  time 
t  =  koT  +  5T  at  the  latest.  If  x(kT)  reaches  the  origin  before  t  =  koT  +  ST,  zero 
inputs  can  be  applied  to  ensure  that  x(kT)  remains  at  the  origin  until  this  time. 
Applying  this  input  sequence  to  the  x-varying  representation. 


z((ko+5)Pr)  = 


h((ko+5)Px) ! 


From  Theorem  3.1,  with  u((ko+5)Px)  =  0, 


z((k0+5)Px+x)  = 


E,  0 
E2  0 


z((k0+5)Px)  =  0. 


Thus,  an  arbitrary  state  zo  has  been  driven  to  0  by  time 

t  =  (no+i+l-t-5P)x  <  nox  +  (5+l)T. 


Theorem  4.3: 


The  x-varying  representation  is  reconstructible  if  the  corresponding 
T-expanded  representation  is  reconstructible,  and  the  state  of  the  corresponding 
representation  at  any  time  nox,  z(noX),  may  be  determined  by  examining  the  inputs 


and  outputs  of  the  system  over  a  period  of  time  no  greater  than  (5+l)T  prior  to 
t  =  n0t. 

Proof: 

Let  no  =  koP  +  j,  j  e  {1,...,P}.  Since  the  T-expanded  representation  is 
reconstructible,  the  state  x(koT)  =  x(koPt)  may  be  determined  by  examining  the 
input  U(kT)  and  output  Y(kT)  of  the  T-expanded  representation  at  the  times 
(ko-6)T,...,(ko-l)T.  From  Theorem  3.1, 


z((koP+l)t)  = 


'Ei  0‘ 

x(koPx) 

.  E2  0  _ 

_h(koPx)_ 

+  B(0)  u(koPt)  = 


'Ef 

.E2. 


x(koT)  +  B(0)  u(koPt). 


Thus,  z((koP+l)t)  may  be  determined  from  knowledge  of  x(koT).  It  follows  that 
z((koP+j)t)  =  z(not)  may  be  determined  for  any  j  e  { 1....P}  by  examining  the 
inputs  and  outputs  of  the  system  at  times  (ko-5)Px,...,(koP+i-l)x  (recall  that  Y(kT) 
involves  values  of  y  at  and  after  time  kPr.)  Therefore,  the  x-varying  representation 
is  reconstructible  and  its  present  state  can  be  determined  by  examining  inputs  and 
outputs  of  the  system  over  a  period  of  time  (nox  -  (ko-5)Px)  <  (5+1  )T  prior  to 
t  =  n0x. 

♦  ♦  ♦ 

Besides  the  application  to  determining  the  controllability  or  reconstructibility  of  the 
t-varving  representation  of  a  multirate  system,  Theorems  4.2  and  4.3  contain  a  subtler  result. 
Results  in  Grasselli  (1984)  show  that  a  x-varying,  T-periodic,  discrete-time  linear  system  of  order 
o  +  rj,  such  as  the  x-varying  representation  under  consideration,  can  be  driven  from  z(nox)  to  the 
zero  state  by  at  most  time  t  =  not  +  (5  +  r|)T  if  it  is  controllable  and  that  the  state  of  the  system  at 
time  t  =  not  can  be  determined  from  knowledge  of  its  inputs  and  outputs  over  a  period  of  time  no 
greater  than  (5  +  r|)T  prior  to  t  =  not  if  it  is  reconstructible.  From  the  special  properties  of 
x-varying  representations  and  their  corresponding  T-expanded  representations.  Theorems  4.2  and 
4.3  shorten  these  time  bounds  by  (rj- 1  )T. 

4.3  Duality,  Stabilizability,  and  Detectability 

In  analogy  with  time-invariant  discrete-time  state  space  representations,  the  reader  may 
anticipate  the  existence  of  stabilizability  and  detectability  properties  and  dual  relationships  between 
stabilizability  and  detectability  for  periodic  discrete-time  representations.  By  combining  results 
from  Grasselli  (1984),  Grasselli  and  Lampariello  (1981),  and  Weiss  (1972),  this  is  seen  to  be  the 
case.  In  the  remainder  of  this  section,  the  argument  x  will  be  suppressed. 
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The  approach  taken  in  defining  stabilizability  and  detectability  will  be  through  a 
decomposition  based  on  controllability  and  reconstructibility.  The  following  lemma  is  a  partial 
restatement  of  results  in  Grasselli  (1984). 

Lemma  4.1: 

A  T-periodic  coordinate  transformation  z(n)  =  A(n)z(n),  A(n)  nonsingular 
for  each  n,  exists  such  that  if  £(n)  =  (A(n),  B(n),  C(n),  D(n))  is  the  representation 
of  Z(n)  in  the  new  coordinates, 


then 


£(n)  =  (A(n+l)A(n)A'Hn),  A(n+l)B(n),  C(n)A-l(n),  D(n)), 


”  au(n) 

ai2(n) 

at3(n) 

a14(n)  “ 

bi(n) 

A(n)  = 

0 

a22(n) 

0 

a24(n) 

,  B(n)  = 

b^n) 

0 

0 

a33(n) 

a34(n) 

0 

_  0 

0 

0 

a44(n)  _ 

-  0  - 

C(n)  =  [  0  c2(n)  0  C4(n)  ]  ,  and  B(n)  =  D(n). 


(4.1) 


All  of  the  submatrices  appearing  in  (4.1)  are  T-periodic  and  have  dimensions  which 
are  constant  with  n.  In  addition,  an(n),  a33(n),  and  a44(n)  are  square  and 
nonsingular  for  all  n. 


an(n)  ai2(n) 
0  a22(n) 

is  controllable,  and 


f  r  a22(n)  a24(n) 


bi(n) 

* 

-Mn). 

,  [  0  c2(n)  ]  ,  B(n) 


V 


0  a44(n)  J 


Mn) 

.  0  . 

,  [  c2(n)  c4(n)  ]  ,  B(n) 


is  reconstructible. 


♦  ♦  ♦ 

As  a  consequence  of  the  results  stated  in  Lemma  4.1,  the  subsystems  (an(n),  bi(n),  0,  D(n)). 
1  a22(n),  b2(n),  c2(n),  D(n)),  (  a33(n),  0,  0,  D(n)),  and  (  aualn),  0,  c?4(n),  D(n))  are.  respective¬ 
ly.  controllable  and  unreconstructible,  controllable  and  reconstructible,  uncontrollable  and 
unreconstructible,  and  uncontrollable  and  reconstructible.  By  the  nonsingularity  of  an(n),  T?33t  n ), 
and  a44(n),  any  zero  eigenvalues  of  (A(P-1)...A(1)A(0)]  must  appear  in  the  controllable  and 
reconstructible  subsystem.  In  light  of  Theorem  3.1,  the  controllable  and  reconstructible  portion  of 
the  t-varying  representation  has  dimension  of  at  least  rj.  This  property  may  serve  to  thwart 
attempts  to  lower  the  dimension  of  the  representation  by  discarding  a'l  subsystems  except  the 
controllable  and  reconstructible  subsystem;  the  remaining  system  will  have  dimension  of  at  least  rj. 
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It  is  suspected  that  a  decomposition  based  on  reachability  and  observability  would  avoid  this 
problem;  however,  the  dimensions  of  subsystems  in  such  a  decomposition  in  general  vary  with 
dme. 

Based  on  the  decomposition  in  Lemma  4.1,  stabilizability  and  detectability  can  be  defined. 

Definition  4.1: 

X(n)  will  be  said  to  be  stabilizabie  if  its  uncontrollable  portion  is 
asymptotically  stable;  with  reference  to  the  decomposition  in  Lemma  4.1,  the 
eigenvalues  of  [  233(P-1).. .233(1)  233(0)]  and  of  [  a^P-l). ..244(1)  244(0)]  all 
have  magnitudes  less  than  1. 

♦  ♦  ♦ 

Definition  4.2: 

X(n)  V'U:  be  said  to  be  detectable  if  its  unreconstructible  portion  is 
asymptotics >:/  sr^ble;  with  reference  to  the  decomposition  in  Lemma  4.1,  the 
eigenvalues  of  [  233(P-1)... 233(1)  233(0)]  and  of  [  2n(P-l)...2n(l)  2n(0)]  all 
have  magnitudes  less  than  1. 

♦  ♦  ♦ 

The  following  lemma,  which  arises  from  results  in  Weiss  (1972)  and  the  periodicity  of  Z(n),  was 
noted  in  Grasselli  and  Lampariello  (1981). 

Lemma  4.2: 

X(n)  =  (A(n),B(n),C(n),D(n))  is  controllable  (reconstructible)  if  and  only  if 
the  dual  of  X(n),  Zd(n)  =  (A,(-n),B,(-n),C'(-n),D,(-n)),  is  reconstructible 
(controllable). 

♦  ♦  ♦ 

By  making  use  of  the  decomposition  in  Lemma  4.1  and  the  duality  properties  in  Lemma  4.2, 
stabilizability  and  detectability  can  be  shown  to  be  dual  properties. 

Theorem  4.4: 

X(n)  is  stabilizabie  (detectable)  if  and  only  if  Xd(n)  is  detectable  (stabilizabie.) 

Proof: 

Let  the  dual  representation  have  state  £(n)  and  let 
Zd(n)  =  (A’(-n),  B’(-n),  C'(- n),  D’(-n))  =  (cc(n),  (3(n),  y(n),  5(n)). 

Consider  the  coordinate  transformation  £(n)  =  ©(n)£(n),  where  @(n)  =  (A*l(  1-n)]'. 

The  representation  of  Xd(n)  in  these  coordinates,  £d(n),  is  given  by 

otfn)  =  @(n+l)  a (n)  ©-^n)  =  [A^f-n)]'  A’(-n)  [A(l-n)]'  =  [A(-n)]\ 

P(n)  =  ©(n+1)  (3(n)  -  [A-'(-n)]'  C(-n)  =  [C(-n)]\ 


and  similarly,  y(n)  =  [B(-n)]'  and  S(n)  =  [5  (-n)]'.  Thus, 


*  * 


a(n)  = 


an(-n)  0 

ai2(-n)  a22(-n) 
ai3(-n)  0 


a33(-n) 


ai4(-n)  a24(-n)  a34(-n)  a^f-n) 


,  P(n)  = 


c2(-n) 


c4(-n), 


y(n)  =  [ "Bl (-n)  b2(-n)  0  0],  and  S(n)  =  D'(-n). 

From  the  properties  of  the  submatrices  in  (4.2)  given  in  Lemma  4.1  and  from 
Lemma  4.2,  it  can  be  seen  that  a33('n)  and  a^-n)  determine  the  dynamics  of  the 
unreconstructible  portion  of  £  d(n)  and  a’n(-n)  and  a33(-n)  determine  the 
dynamics  of  the  uncontrollable  portion  of  2d(n).  Note  that 

au(-(P-l))...an(-l)  an(O)  =  an(l)...2n(P-l)  an(0) 

=  [an(0)a11(P-l)..a11(l)l'. 

Sincea'n(.)  is  square,  the  eigenvalues  of  a'n(-(P-l))...an(- 1  )an(0)  are  equal 
to  the  eigenvalues  of  an(P-l)...au(l)  an(0).  A  similar  relation  holds  fora33(-n) 
and  ^(-n).  Thus,  the  uncontrollable  (unreconstructible)  portion  of  £(n)  is 
asymptotically  stable  if  and  only  if  the  unreconstructible  (uncontrollable)  portion  of 
2d(n)  is  asymptotically  stable.  The  statement  of  the  theorem  then  follows  from 
Definitions  4.1  and  4.2. 

♦  ♦  ♦ 

The  examination  of  the  relationship  between  corresponding  T-expanded  and  x-varying 
representations  can  now  be  completed. 

Lemma  4.3: 

(A(n),  B(n),  C(n),  D(n))  is  stabilizable  if  there  exist  T-periodic  matrices 
F(n)  such  that  on  applying  the  state  feedback 

u(kP+i)  =  -F(kP+i)  z(kP),  0  <  i  <  P,  k  =  0,1,...  , 

the  resulting  representation  is  asymptotically  stable;  that  is,  the  closed-loop 
monodromy  matrix  at  time  t  =  0, 

<PC(0)  =  A(P-l)  [,..[A(1)  [A(0)  -  B(0)  F(0)]  -  B(l)  F(l)]„.]  -  B(P-l)  F(P-l), 


Proof: 


has  eigenvalues  of  magnitude  less  than  1. 

Without  loss  of  generality,  it  can  be  assumed  that  (A(n),  B(n),  C(n),  D(n)) 
has  been  transformed  into  the  canonical  form  in  Lemma  4. 1 ;  it  can  be  verified  that  if 


state  feedback  P(n)  =  F(n)  A(0)  is  used  in  the  z  coordinates,  then 

<$c(0)  =  A'HP)  <I>c(0)  A(0)  =  A-HO)  <DC(0)  A(0). 

Thus,  $c(0)  and  <t>c(0)  have  the  same  eigenvalues.  Partition  (A(n),  B(n))  into  its 
controllable  and  uncontrollable  pans: 


A(n)  = 


An(n)  A12(n) 

0  A22(n) . 


Note  that  if  there  is  F(n)  such  that  each  eigenvalue  of 

<$c(0)  =  A(P-i)  [...[A(i)[A(0)-B(0)F(0)]-B(i)F(i)]...]  -  B(P-i)F(P-i) 

has  magnitude  less  than  1,  then  each  eigenvalue  of  [A22(P-1)...A22(1)  A22(0)]  has 
magnitude  less  than  1.  This  is  easily  seen  by  expanding  <5C(0)  as 

P-2 

ac(0)  =  A(P-1)...A(1)  A(0)  -  z  [A(P-l)...A(j-t-l)S(j)Ti<j)]  -  B(P-1)  F(P-l) 

j=0 


and  observing  that  A(P-l)...A(j+l)  is  upper  triangular  for  -1  <  j  <  P-2.  From  the 
special  form  of  B(n),  <J>C(0)  is  upper  triangular  with  [A22(P-1)...A22(1)  A22(0)] 
appearing  on  the  diagonal.  Thus,  if  the  eigenvalues  of  <5C(0)  have  magnitudes  less 
than  1,  so  do  the  eigenvalues  of  [A22(P-1)...£22(1)  A22(0)],  and  (A(n),  B(n))  is 
stabilizable. 

♦  ♦  ♦ 

Theorem  4,5: 

The  x-varying  representation  of  a  multirate  system  is  stabilizable  if  its 
corresponding  T-expanded  representation  is  stabilizable. 

Proof: 

The  theorem  can  be  restated  equivalently  as  “the  x-varying  representation  of 
a  multirate  system  is  stabilizable  if  its  completely  corresponding  T-expanded 
representation  is  stabilizable,”  since  the  input  matrix  of  the  completely  correspond¬ 
ing  representation,  Bcc,  is  merely  a  rearrangement  of  the  columns  of  the  input 
matrix  of  the  corresponding  representadon,  Be,  padded  with  columns  of  zeroes. 

If  the  completely  corresponding  T-expanded  representation  is  stabilizable,  it 
is  a  standard  time-invariant  result  that  there  is  Fcc  such  that  the  eigenvalues  of 
(Ac  -  BCCFCC)  have  magnitudes  less  than  1.  Consider  applying  state  feedback  to 
the  x-varying  representation  as  in  Lemma  4.3: 


u(kP+i)  =  -F(kP+i)  z(kP),  0  <  i  <  P,  k  =  0,1 . 

Suppose  u  €  aPx1.  Partition  Fcc  into  P  groups  of  rows. 
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i-  , 


Fcc  = 


where  f;  e  RPx5.  Set  F(n)  =  [  f„  Opx^]  and  write  Oc(0)  as 


<DC(0)  =  A(P-1)...A(1)  A(O)  -  X  [A(P-l)...A(j+l)B(j)F(j)]  -  B(P-l)  F(P-l). 

j=0 


In  the  notation  of  Theorem  3.2, 


<DC(0)  =  <x  -  (3 


F(P-l) 


From  the  values  of  F(n)  and  Theorem  3.2, 


r  Ae  0  1  [Bcc" 

^c(O)  =  Q  '  J  ^  ^CC 


Thus,  the  eigenvalues  of  Oc(0)  have  magnitudes  less  than  1,  and  by  Lemma  4.3, 
the  T-varying  representation  is  stabilizable. 


By  using  Theorem  3.2,  Theorem  4.4,  and  Lemma  4.3  and  proceeding  in  a  manner  parallel  to 
Theorem  4.5,  the  following  theorem  can  also  be  established. 

Theorem  4.6: 

The  T-varying  representation  of  a  multirate  system  is  detectable  if  its 
corresponding  T-expanded  representation  is  detectable. 


The  statements  of  Theorems  4.2, 4.3, 4.5,  and  4.6  all  hold  with  “if’  replaced  by  “only  if,” 
but  the  new  assertion  is  rather  unimportant  and  in  general  trivial  to  prove.  By  arguments  parallel  to 
those  presented  in  this  chapter,  it  can  be  shown  that  M-varying  representations  which  correspond 
to  T-expanded  representations  (in  a  manner  analogous  to  corresponding  T-varying  and  T-expanded 
representations)  inherit  the  stability,  controllability,  reconstructibility,  stabilizability,  and  detect¬ 
ability  properties  of  the  T-expanded  representation. 


'V.VLV 


CHAPTER  5 


DESIGN  OF  MULTIRATE  CONTROLLERS 


The  representations  and  methods  of  analysis  developed  for  multirate  systems  thus  far  bear 
fruit  in  the  form  of  multirate  controllers.  The  term  “multirate  controller”  need  not  imply  that  only 
the  controller  is  multirate;  candidates  for  control  include  all  systems  to  which  Procedure  2.4  and  its 
extensions  apply.  Although  the  T-expanded  or  M-varying  representations  are  not  truly  time 
invariant  or  periodic  due  to  the  time  expanded  form  of  their  inputs  and  outputs,  this  characteristic 
offers  little  or  no  obstruction  to  the  design  and  implementation  of  multirate  controllers  by 
conventional  means.  The  controller  designs  presented  here  are  intended  as  illustrations  of  the 
general  properties  of  controller  design  and  implementation  using  T-expanded  and  M-varying 
representations  of  multirate  systems.  The  usefulness  of  these  representations  is  by  no  means 
limited  to  the  controllers  examined  here. 


5.1  T-expanded  State  Feedback  and  Observers 

Consider  a  multirate  system  with  the  T-expanded  representation 
x((k+l)T)  =  Ae  x(kT)  +  Bve  V(kT)  +  Bre  R(kT) 

Y(kT)  =  Ce  x(kT)  +  Dve  V(kT)  +  Dre  R(kT),  ' " ' '  * 

where  V(kT)  is  an  input  intended  for  control  purposes  and  R(kT)  is  a  reference  or  load  input.  On 
applying  feedback  of  the  form  V(kT)  =  -F  x(kT),  the  state  equation  becomes 

x((k+l)T)  =  [Ae  -  BveF]  x(kT)  +  BreR(kT). 

If  a  matrix  F  can  be  found  such  that  the  eigenvalues  of  [Ae  -  BveF]  are  all  less  than  1  in  magnitude, 
then  Theorem  4.1  asserts  that  this  representation,  as  well  as  the  corresponding  t-varying 
representation,  is  asymptotically  stable.  Conditions  for  the  existence  of  such  a  matrix  F  given 
(Ae,  Bve)  are  well  known,  and  methods  of  calculating  F  abound. 

Although  calculation  of  F  by  treating  (Ae,  Bve)  as  a  time- invariant  pair  using  pole 
placement  or  linear  quadratic  methods  serves  the  purpose  of  stabilizing  the  system,  two  comments 
are  in  order  regarding  the  choice  of  design  parameters  for  these  methods.  When  choosing  pole 
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locations  for  pole  placement,  bear  in  mind  that  these  pole  locations  are  in  a  time  scale  of  T  seconds 
per  transition.  An  LQ  design  with  even  weighting  on  each  member  of  an  expanded  bundle  of 
inputs  with  period  Ti,  say  Vi(kT),  may  result  in  a  control  law  which  consistently  produces  inputs 
v^kT)  of  magnitude  quite  different  from  vi((k+l)T-T| ).  The  cause  of  this  behavior  is  that  if  the 
subsystems  affected  by  vi  have  poles  far  from  the  origin  (close  to  the  origin,)  then  an  input  V](kT) 
will  have  a  greater  (lesser)  effect  on  x((k+l)T)  than  an  input  vi((k+l)T-Ti).  If  such  behavior  is 
deemed  undesirable,  it  can  be  predicted  by  examining  the  columns  of  Bve  corresponding  to  V  i(kT) 
and  corrected  by  appropriate  changes  in  the  input  weighting  matrix. 

The  design  of  observers  which  produce  an  estimate  x(kT)  of  x(kT)  parallels  the  design  of 
such  observers  in  the  standard  time-invariant  case,  with  two  minor  exceptions.  The  observer  must 
produce  a  predictive  estimate,  usually  denoted  by  x(kT  I  (k-l)T),  since  current  estimates, 
commonly  denoted  by  x(kT  I  kT),  result  in  a  noncausal  observer.  The  observer  must  also  account 
for  the  Dve  and  Dre  terms  frequently  present  in  the  T-expanded  representation. 

The  state  equation  for  an  observer  for  system  (5.1)  in  predictive  form  is 

x((k+l)T)  =  Aex(kT)  +  BveV(kT)  +  BreR(kT)  +  K[Y(kT)  -  Cex(kT)  -  DveV(kT)  -  DreR(kT)]. 

The  estimation  error  satisfies  x((k+l)T)  -  x((k+l)T)  =  [Ae  -  KCe](x(kT)  -  x(kT)).  If  K  is  chosen  to 
stabilize  [Ae  -  KCe]  by  methods  dual  to  those  employed  to  find  F,  this  error  will  asymptotically 
approach  zero  as  k— ►«>.  The  state  and  output  equations  for  the  combined  state  feedback  and 
observer  pair  are 

x((k+l)T)  =  [Ae  -  KCe  *  [Bve  -  KDve)F]  x(kT)  +  KY(kT)  +  [Bre  -  KDre]R(kT) 

(5.2) 

V(kT)  =  -F  x(kT). 

It  is  easily  verified  that  the  principle  of  separation  holds  for  the  state  feedback  and  observer 
designs.  Equation  (5.2)  is  a  multirate  discrete-time  system  specified  in  T-expanded  form  with 
input  [Y(kT)'  R(kT)']'  and  output  V(kT).  As  x((k+l)T)  is  needed  at  time  t  =  (k+l)T  to  compute 
V((k+1)T)  and  each  component  of  Y(kT)  and  R(kT)  is  available  at  time  t  =  ((k+l)T  - 1)  at  the 
latest,  it  follows  that  the  controller  given  above  is  causal.  Although  all  of  the  information 
necessary  to  calculate  x((k+l)T)  may  not  be  present  until  time  t  =  ((k+l)T-  x),  the  calculation  of 

[Ae  -  KCe  -  [Bve  -  KDve]F]  x(kT) 
is  possible  immediately  after  t  =  kT.  From 

K  Y  (kT)  =  Kuyi(kT)  +  Kii2yi(kT+Ti)  +...+  K2>1y2(kT)  +  K2<2y2(kT-rT2)  +..., 

where  the  Kj,j  are  groups  of  columns  of  K  chosen  to  be  compatible  with  the  individual  values  of 
the  i!h  bundle  of  signals  (with  period  T;)  comprising  Y(kT)  at  time  t  =  kT  +  jTi,  the  calculation  of 
KY(kT)  and  [Bre  -  KDre|R(kT)  can  occur  progressively  as  the  values  of  Y(kT)  and  R(kT)  become 
available.  The  implementation  of  (5.2)  involves  sampling  each  component  of  v(t)  and  r(t)  at  the 


appropriate  rate  and  transferring  values  of  -F  x(kT)  to  each  component  of  v(t  ;  at  the  appropriate 
rate  and  holding  each  of  these  values  until  the  next  arrives.  This  behavior  is  conceptualized  as 
sample/hold  devices  at  appropriate  rates  at  the  inputs  and  outputs  of  the  controller.  The  following 
procedure  summarizes  the  steps  necessary  to  find  a  T-expanded  multirate  state  feedback- 
observer  controller. 

Procedure  5.1: 

a.  Decide  upon  the  rate  desired  for  each  of  the  controller  outputs  and  the 
desired  rate  at  which  each  plant  output  and  reference  input  is  to  be  observed  by  the 
controller.  Place  samplers  of  corresponding  periods  on  the  block  diagram  of  the 
system  after  the  controller  output  v  and  the  reference  input  r  and  before  the 
controller  input  y,  as  shown  in  Figure  5.1.  Assume  that  only  the  components  of  r 
observed  by  the  controller  are  fed  to  the  summer,  as  in  Figure  5.1.  Relaxation  of 
this  assumption  is  considered  after  this  procedure. 


Figure  5. 1 .  Placement  of  samplers. 


b.  Apply  Procedure  2.4  and  any  modifications  which  apply  to  find  a 
T-expanded  representation  such  as  (5.1),  treating  r  and  v  as  inputs  and  y  as  the 
output. 

c.  Using  any  method  desired,  find  matrices  F  and  K  from  the  pairs 
( Ae,  Bve)  and  (Ae,  Ce)  such  that  each  eigenvalue  of  [Ae  -  BveF]  and  [Ae  -  KCc] 
has  magnitude  less  than  1.  The  controller  state  and  output  equations  are  given  by 

(5.2). 

♦  ♦  ♦ 

The  assumption  in  Procedure  5.1(a)  that  only  the  components  of  the  reference  input 
observed  by  the  controller  are  fed  to  the  summer  is  necessary  for  the  design  procedure  but  need  not 
be  satisfied  in  the  implementation  of  the  controller.  Let  rs  be  the  actual  sampled  and  held  version  of 
r  applied  to  the  plant,  where  each  component  of  rs  has  a  period  which  is  an  integer  multiple  of  t. 
the  fundamental  period  used  in  the  execution  of  Procedure  5.1.  Let  r0  be  the  sampled  version  of  r 
observed  by  the  controller.  Since  the  components  of  rs  and  r0  vary  at  times  which  are  integer 
multiples  of  t,  rs  can  be  written  rs  =  rQ  +  ru,  where  the  components  of  ru  each  vary  with  period  t. 


The  signal  ru  is  the  portion  of  rs  unobserved  by  the  controller.  Let  the  plant  and  controller  system 
with  r0  and  ru  as  inputs  have  composite  state 

r  x(kT)  *] 


and  state  equation 


s(kT)  = 

Lx(kT). 

s((k+l)T)  =  A  s(kT)  +  BoRo(kT)  +  BuRu(kT). 


(5.3) 


A  is  stable  by  design;  therefore,  under  the  quite  reasonable  assumption  that  B0  and  Bu  are 
bounded,  (5.3)  is  bounded-input  bounded-state  stable.  If  it  is  desired  that  the  controller  sample  r  at 
a  slow  rate,  then  either  r  must  be  fed  to  the  plant  at  this  slow  rate,  slowing  the  system's  response 
time,  or  r  can  be  fed  to  the  plant  at  a  fast  rate,  in  which  case  the  unobserved  portion  ru  acts  as  a 
disturbance 

The  T-expanded  multirate  controller  design  has  the  advantages  of  ease  of  calculation  by 
established  methods  and  being  readily  implemented.  To  its  disadvantage,  this  controller  exhibits  a 
ratiier  sluggish  response  to  rapid  changes  in  r  if  T  is  large.  The  state  feedback  V(kT)  =  -F  x(kT)  is 
essentially  applied  open  loop  over  times  of  length  T,  leaving  the  system  vulnerable  to  disturbances 
and  inputs  ru  unobserved  by  the  controller.  Indeed,  if  the  observer  portion  is  initialized  to  x(0)  =  0. 
no  control  will  be  applied  until  time  t  =  T.  The  design  of  a  T-expanded  multirate  Kalman  filter  for 
a  system  with  noise  at  the  input  which  is  uncorrelated  with  the  noise  at  the  output  is  complicated  by 
the  fact  that  the  De  term  in  the  T-expanded  representation  introduces  noise  at  the  output  which  is 
highly  correlated  with  the  noise  at  the  input. 


5.2  M-varying  State  Feedback  and  Observers 

With  the  intent  of  achieving  a  faster  response  to  inputs  and  disturbances,  the  value  of  the 
state  may  be  estimated  and  fed  back  at  times  t  =  nx.  Consider  the  x-varying  representation  of  the 
multirate  system  in  Figure  5.1  (in  the  following,  the  parameter  x  is  suppressed): 

z(n+l)  =  A(n)  z(n)  +  Bv(n)  v(n)  +  Br(n)  r(n) 

(5.4) 

y(n)  =  C(n)  z(n)  +  Dv(n)  v(n)  +  Dr(n)  r(n). 

After  applying  state  feedback  v(n)  =  -F(n)  z(n), 

z(n-t-l)  =  (A(n)  -  Bv(n)F(n)]  z(n)  +  Br(n)  r(n). 

If  F(0),...,  F(P-l)  can  be  found  such  that  each  eigenvalue  of  the  closed-loop  monodromy  matrix  at 
t  =  0, 

<t\(0)  =  [A(P-l)  -  BV(P-1  )F(P- 1 )]...[ A(  1 )  -  Bv(  1)F(  1 )  1  [ A( 0)  -  Bw(0)F(0)]. 
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has  magnitude  less  than  1,  then  the  system  employing  F(n)  =  F(n  mod  P)  as  T-penodic  feedback 
matrices  will  be  asymptotically  stable.  Such  a  periodic  feedback  arises  naturally  as  the  solution  to 

the  LQ  problem  with  the  cost  function 


J  =  1/2  ^  [z'(n)Wz'(n)Wz(n)z(n)  +  v'(n)Wv'(n)Wv(n iv< n ij,  (5.5) 

n=0 

where  Wz(n+P)  =  Wz(n)  and  Wv(n+P)  =  Wv(n)  for  each  n  e  Z+.  The  discrete-time  Riccati 
equation  associated  with  this  problem  admits  aT-periodic  solution  which  determines  the  feedback 
matrices  F(n).  Discrete-time  periodic  Riccati  equations  have  just  recently  been  investigated; 
conditions  for  the  existence  and  uniqueness  of  their  solutions,  methods  of  calculating  these 
solutions,  and  conditions  under  which  the  resulting  feedback  is  stabilizing  are  still  developing.  See 
Bittani  et  al.  (1986)  and  the  references  cited  therein. 

In  order  to  make  use  of  the  results  in  Bittani  et  al.  (1986),  assume  that  Wv(n)  is  square  and 
nonsinguiar  for  each  n  €  Z+.  The  Riccati  equation  associated  with  (5.4)  and  (5.5)  can  then  be 
written 

S(n)  =  A'(n)S(n+l)A(n)  +  Wz'(n)Wz(n) 

-  A'(n)S(n+l)Bw(n)  [I  +  Bw’(n)S(n+l)Bw(n)]-1  Bw’(n)S(n+l  )A(n),  (5.6) 

where  Bw(n)  =  Bv(n)(Wv(n))'1.  For  state  feedback  v(n)  =  -F(n)  z(n), 

F(n)  =  (Wv(n))-1  [I  +  Bw'(n)S(n+l)Bw(n))-1  Bw’(n)S(n+l)A(n) 

=  [Wv’(n)Wv(n)  +  Bv'(n)S(n+l)Bv(n)]-1  Bv'(n)S(n+l)A(n).  (5.7) 

The  following  theorem  is  stated  without  proof;  its  proof  relies  on  the  application  of  results  in 
Bittani  et  al.  (1986),  the  duality  properties  of  discrete-time  periodic  systems  (see  Theorem  4.4)  and 
the  assumed  nonsingularity  of  Wv(n). 

Theorem  5.1: 

Consider  the  dynamical  system  represented  by 

z(n+l)  =  A(n)  z(n)  +  Bv(n)  v(n) 

(5.8) 

y(n)  =  Wz(n)  z(n), 

where  A(n)  and  Bv(n)  are  as  in  (5.4)  and  Wz(n)  is  as  in  (5.5).  Then  a  unique 
T-periodic  symmetric  positive  semidefinite  solution  to  (5.6)  exists  and  (5.4)  is 
asymptotically  stable  after  applying  feedback  v(n)  =  -F(n)  z(n),  where  F(n)  is  given 
by  (5.7),  if  and  only  if  (5.8)  is  stabilizable  and  detectable;  i.e.,  the  uncontrollable 
part  of  (5.8)  is  asymptotically  stable  and  the  unreconstructible  pan  of  (5.8)  is 
asymptotically  stable. 

♦  ♦  ♦ 

By  Theorem  4.5,  the  stabilizability  of  (5.8)  can  be  evaluated  by  examining  the  stabilizability  of  the 


T-expanded  representation  (with  r(n)  =  0)  corresponding  to  (5.4).  This  device  cannot  be  employed 
to  determine  whether  (5.8)  is  detectable;  however,  if  Wz(n)  is  square  and  nonsingular  for  each  n, 
(5.8)  is  detectable. 

A  t-varving  observer  to  form  an  estimate  z(n)  of  z(n)  takes  a  form  identical  to  the 
T-expanded  observer: 

z(n+l)  =  A(n)z(n)  +  Bv(n)v(n)  +Br(n)r(n)  +  K(n)[y(n)  -  C(n)z(n)  -  Dv(n)v(n)  -  Dr(n)r(n)]. 

Since  z(n+l)  -  z(n+l)  =  [A(n)  -  K(n)C(n)]  (z(n)  -  z(n)),  the  error  z(n)  -  z(n)  will  asymptotically 
approach  zero  if  K(n)  is  chosen  in  a  manner  dual  to  the  choice  of  F(n)  for  stabilization  (see 
Theorem  4.4  for  the  required  duality  details.)  By  writing  the  equations  of  the  combined  controller 
and  plant  in  terms  of  z(n)  and  (z(n)  -  z(n))  and  Finding  the  T-expanded  representation  of  this 
T-varying  representation,  it  can  be  verified  that  the  principle  of  separation  holds  for  the  state 
feedback  and  observer  designs. 

The  procedure  for  finding  a  T-varying,  T-periodic  state  feedback-observer  controller  via  the 
discrete-time  periodic  Riccati  equation  parallels  Procedure  5.1.  Comments  similar  to  those 
following  Procedure  5.1  with  regard  to  reference  inputs  unobserved  by  the  controller  also  apply. 

Using  the  periodic  description  (5.4)  for  the  multirate  system,  other  types  of  T-varying 
controllers  can  be  designed.  Let  z(n)  €  R5+Tl,  as  in  Chapter  4.  Results  in  Grasselli  and 
Lampariello  (1981)  show  that  the  reconstructibility  and  controllability  of  (5.4)  (with  r(n)  =  0)  are 
necessary  and  sufficient  for  the  existence  of  a  T-varying,  T-periodic  state  feedback-observer  dead¬ 
beat  controller  for  (5.4)  which  is  capable  of  driving  any  initial  condition  to  the  origin  within  an 
interval  of  time  no  greater  than  2(8+r|)T.  In  light  of  Theorems  4.2  and  4.3,  it  is  suspected  that  a 
closer  examination  of  the  dead-beat  controller  problem  could  reduce  this  interval  to  2(5+1  )T  for  the 
representation  (5.4).  M-varying  controllers  can  be  found  by  utilizing  the  M-varying  representation 
of  the  multirate  system  in  Figure  5.1.  As  M-varying  representations  involve  M-expanded  inputs 
and  outputs,  care  must  be  taken  to  insure  that  the  controller  designed  is  causal.  If  the  M-varying 
controller  is  of  the  state  feedback-observer  type,  the  controller  will  be  causal  if  the  observer 
produces  predictive  estimates. 

Implementation  of  a  T-varying  controller  requires  an  interpretation  of  v(n),  r(n),  and  y(n). 
In  Procedure  5.1(a),  sample/hold  devices  are  placed  on  the  block  diagram  to  represent  the  effects 
of  the  sampling  and  “output  and  hold”  operations  performed  by  the  discrete  controller.  The  fact 
that  values  of  v,  r,  and  y  are  specified  for  the  controller  at  times  t  =  me  might  suggest  that  the 
T-varying  controller  be  implemented  by  loading  values  of  r,  y,  and  v  into  and  out  of  the  controller 
at  times  t  =  nT  and  physically  passing  these  values  through  the  sample/hold  devices  following  v 
and  r  and  preceding  y  on  the  block  diagram  of  Procedure  5.1(a).  The  implementation  can  be  made 
much  more  efficient  by  bringing  the  sample/hold  devices  inside  the  controller  in  the  following 
manner.  For  each  component  Vj(n)  of  v(n).  if  n  corresponds  to  a  time  at  which  the  sample/hold 


associated  with  v;  samples,  then  Vj(n)  should  be  output  by  the  controller  at  time  t  =  nx  and  held  for 
the  duration  of  the  sampling  period  associated  with  v^  For  each  component  Vj(n)  of  y(n),  if  n 
corresponds  to  a  time  at  which  the  sample/hold  associated  with  yj  samples,  then  yj(n)  should  be 
acquired  by  the  controller  at  time  t  =  nx  and  stored  for  the  duration  of  the  sampling  period 
associated  with  yj  to  provide  the  values  yj(n),  yj(n+l),... .  The  signals  r(n)  are  treated  similar  to 
v(n).  Analogous  comments  hold  for  M-varying  controllers. 

A  x-varying  controller  can  potentially  respond  quickly  to  changes  in  the  input  r  and,  since 
the  state  estimate  is  updated  frequently,  reduce  the  effects  of  unmeasured  disturbances  on  the 
output.  The  price  paid  for  these  attributes  takes  the  form  of  complexity  of  design  and 
implementation;  the  x-varying  representation  employed  by  the  controller  has  more  states  than  the 
T-expanded  representation,  resulting  in  greater  computation  and  storage  requirements.  An 
M-varying  controller,  where  M  is  chosen  such  that  T/M  e  N,  may  serve  as  an  adequate 
compromise  between  x-varying  and  T-expanded  controllers;  by  choosing  different  values  of  M, 
complexity  can  be  traded  for  response  time.  The  practicality  of  solving  the  discrete-time  periodic 
Riccati  equation  or  Finding  the  feedback  matrices  required  in  Grasselli  and  Lampariello  (1981)  is 
unknown  as  of  the  present. 

5.3  Transfer  Function  Controller  Designs 


$ 
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By  thinking  of  (5. 1)  as  a  time-invariant  state  and  output  equation,  a  transfer  function  matrix 
may  be  obtained  for  this  system: 

rV(z)1  rV(z)l 

Y(z)  =  (  Ce(zl  -  A e)-1  [Bve  Brc]  +  [Dve  Dre]  )  N  =  G(z)  . 

L  R(z)  J  LR(z)J 

Using  MIMO  transfer  function  design  techniques,  controllers  of  a  variety  of  structures  may  be 
designed  (Kailath,  1980).  One  such  structure  is 

V(z)  =  Vy(z)  +  Vr(z)  =  Hy(z)  Y(z)  +  Hr(z)  R(z). 

Although  details  of  such  a  design  will  not  be  elaborated  upon  here,  the  restrictions  imposed  on 
Hy(z)  and  Hr(z)  by  the  T-expanded  nature  of  V(kT),  Y(kT),  and  R(kT)  warrant  discussion.  A 
state  space  realization  of  Hy(z),  if  it  exists,  takes  the  form  of  a  multirate  discrete-time  system 
specified  in  T-expanded  form, 

w((k+l)T)  =  Ay  w(kT)  +  By  Y(kT) 

(5.9) 

Vy(kT)  =  Cy  w(kT)  +  Dy  Y(kT), 

where  Dy  =  Hy(<»).  For  this  system  to  be  physically  realizeable,  the  elements  of  Hv(°°)  must  be 
Finite;  Hy(z)  must  be  proper.  For  (5.9)  to  be  causal,  selected  entries  in  Dy  must  be  zero,  and  the 


corresponding  entries  in  Hy(z)  must  be  strictly  proper.  Let  Hy(z)  be  partitioned  into  i  x  j  blocks 
hij(z)  to  conform  with  the  signal  bundle  values  v(kT+tj)  comprising  V(kT)  and  y(kT+tj) 
comprising  Y(kT)  (this  notation  is  expedient  and  unrelated  to  previous  notations.)  For  Hv(z)  to  be 
realizable  as  a  causal  multirate  system,  it  must  satisfy 

f  0  ,  if  t,  <  t. 

hij(°°)  =  \  .  .  .  .  (5  10 

(  a  matrix  of  finite  numbers  ,  if  tj  <  tj 

Multirate  controllers  designed  by  transfer  function  tecnniques  possess  some  of  the 
advantages  of  both  the  T-varying  and  t-varying  state  space  controllers.  The  T-expanded  nature  of 
such  a  controller  simplifies  its  implementation,  and  the  possibility  of  incorporating  “D”  terms  into 
the  state  space  realization  of  the  controller  may  allow  the  controller  output  during  [kT,  (k+l)T)  to 
depend  on  its  inputs  during  [kT,  (k+l)T),  avoiding  the  open  loop  type  behavior  of  the  T-expanded 
state  space  controller  with  respect  to  measured  disturbances.  However,  the  state  of  a  multirate 
controller  designed  by  transfer  function  techniques  is  updated  only  at  times  t  =  kT,  which  may 
impair  its  ability  to  reduce  the  effects  of  unmeasured  disturbances  on  the  output.  The  practical 
implications  of  incorporating  the  causality  constraints  (5.10)  into  MIMO  transfer  function  design 
techniques  are  unknown  at  the  present. 

5.4  Comments 

The  advantages  of  using  a  multirate  controller  as  opposed  to  a  single-rate  controller  depend 
greatly  on  the  specific  system  to  be  controlled  and  on  an  intelligent  choice  of  sampling  periods  for 
the  controller.  Clearly,  if  the  plant  is  inherently  multirate,  conventional  methods  of  analysis  and 
design  are  not  applicable  to  controller  design.  However,  multirate  systems  frequently  arise  from 
attempts  to  circumvent  the  “a  controller  is  only  as  fast  as  its  slowest  actuator  or  sensor”  principle  of 
single-rate  sampled-data  designs.  In  such  a  situation,  the  simplicity  of  a  single-rate  controller  must 
be  weighed  against  the  prospects  of  making  full  use  of  the  bandwidth  of  the  actuators  and  sensors 
by  multirate  control.  Results  in  Barnes  and  Shinnaka  (1980)  indicate  that  implementations  of 
multirate  systems  specified  in  T-expanded  form  exhibit  desirable  numerical  characteristics,  such  as 
low  roundoff  noise.  A  basic  property  of  multirate  systems  is  that  as  the  numbers  pi  obtained  from 
the  normalization  process  increase,  representation  and  analysis  of  the  system  become  more 
difficult.  Exploiting  any  freedom  in  the  choice  of  actuator  and  sensor  rates  to  reduce  the  numbers 
Pi  will  result  in  a  much  simpler  multirate  controller. 


CHAPTER  6 


THE  VARIABLE  COMPONENT  METHOD 
APPLIED  TO  MULTIRATE  DISCRETE-TIME  SYSTEMS 

The  variable  component  method  and  the  method  of  sensitivity  points  have  been  applied  to 
synthesize  and  tune  linear,  time-invariant  controllers  for  linear,  time-invariant  plants;  see  Frank 
(1978),  Kokotovic  (1964),  Kokotovic  (1965),  Hung  (1985),  and  the  references  cited  therein.  In 
this  chapter,  the  variable  component  method  will  be  extended  to  the  class  of  multirate  discrete-time 
systems  to  which  Procedure  2.4  and  its  extensions  apply.  A  general  time-invariant  variable 
component  result  will  be  derived  and  then  related  to  the  T-expanded  representation  of  these 
multirate  systems. 

Consider  a  time-invariant  discrete-time  single-rate  system  with  an  embedded  scalar 
parameter  k  and  scalar  output  y(u(n),k,n)  at  time  n  for  a  given  scalar  input  u(n).  The  variable 
component  method  and  the  method  of  sensitivity  points  are  each  based  on  a  block  diagram 
representation  of  the  system  and  provide  a  means  of  determining  the  output  sensitivity 
function,  9[y(u(n),k,n)]/3k,  by  simulation  or  implementation.  These  two  methods  basically 
proceed  by  injecting  the  input  u(n)  into  the  system  or  a  suitable  model  of  the  system  and  extracting 
selected  signals,  as  determined  by  the  block  diagram,  from  the  system  or  model  and  injecting  these 
signals  (or  possibly  Filtered  versions  of  these  signals)  at  selected  points  of  another  model  of  the 
system.  Alternately,  the  extracted  signals  may  be  stored  and  injected  into  the  actual  system  at  a 
later  time.  In  either  case,  the  response  at  a  selected  point  of  the  two  connected  systems  is  the 
output  sensitivity  function.  In  addition  to  the  actual  system,  the  variable  component  method 
requires  as  many  models  of  the  system  as  there  are  parameters  to  be  independently  varied  if  the 
output  sensitivity  function  for  each  parameter  is  to  be  obtained  simultaneously.  The  method  of 
sensitivity  points  allows  the  output  sensitivity  functions  for  any  number  of  parameters  to  be 
obtained  simultaneously  with  only  one  system-model  pair.  However,  the  method  of  sensitivity 
points  cannot  be  applied  to  MIMO  systems  (Hung,  1985). 

Once  the  output  sensitivity  function  has  been  obtained,  the  response  of  the  system  with 
parameter  value  k  +  Ak,  lAkl  «  Ikl,  and  input  u(n)  can  be  approximated  by 

y(u(n),k+Ak,n)  =  y(u(n),k,n)  +  Ak(d[y(u(n),k,n)]/5k). 


v  nmnm 


The  output  sensitivity  function  thus  provides  information  useful  for  iteratively  adjusting  the 
parameter  k  so  that  the  output  of  the  system  meets  or  approaches  certain  time-domain  criteria  for  a 
given  input  u(n).  Such  adjustments  are  frequently  determined  through  the  solution  of  an 
opdmization  problem  by  minimizing  a  measure  of  the  error  between  the  actual  output  of  the  system 
and  the  desired  output  of  the  system  (see  Kokotovic  (1965)  and  Hung  (1985).)  An  example  of  the 
use  of  the  variable  component  method  with  a  multirate  system  would  be  the  tuning  of  parameters  in 
a  single-rate  sampled-data  controller  based  on  intersample  values  of  the  output  of  the  controlled 
system. 


6.1  A  Time-invariant  Result 

In  the  following,  a  general  variable  component  method  result  for  time-invariant  discrete¬ 
time  MIMO  systems  will  be  derived.  This  time-invariant  result  can  then  be  applied  to  the 
T-expanded  representation  of  a  multirate  discrete-time  system. 

Consider  a  MIMO  time-invariant  single-rate  discrete-time  linear  system  with  rate  1/T  and  a 
parameter  K  =  diag[ki,k2,...,kr]  which  can  be  isolated  as  shown  in  Figure  6.1.  Let  U(z),  E(z), 
and  Y(z)  be  the  vector  Z-transforms  of  U(mT),  E(mT),  and  Y(mT),  respectively,  and  let  U(z)  be 
j  x  1,  E(z)  be  r  x  1,  and  Y(z)  be  q  x  1.  The  dependence  of  E(z)  and  Y(z)  on  K  and  U(z)  will 
not  be  made  explicit  at  this  point. 


Figure  6.1.  A  MEMO  time-invariant  system. 


Assuming  that  the  system  in  Figure  6.1  is  well-posed  (see  Chen  (1984),)  there  exist  transfer 

function  matrices  F(z),  G(z),  H(z),  and  J(z)  which  do  not  depend  on  the  parameters  kj,  k: . kr 

such  that,  suppressing  the  argument  z, 


Thus, 

Also, 


\  *.  *.  *. ’ 
*Culj£  m 


R  =  GU  +  HKE  =  GU  +  HKR. 
R  =  (Ir-  HK)1  GU. 


Y  =  FU  +  JKR  =  FU  +  JK  (Ir  -  HK)1  GU. 
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(6.1 ) 
(6.2) 


For  the  set  S  £  { l,2,...,r},  define  the  r  x  r  matrix  A(r,S)  by 

A  ,  c,  f  1  ,  if  i  =  j  e  S 

Aij(r,S)  =  \ 

l  0  ,  otherwise 

Note  that  dK/dk,  =  A(r,  { i } ).  From  (6.2),  we  can  compute 

9Y/9kj  =  J  A(r, ( i } )  (Ir  -  HK)-1  GU  +  JK  (Ir  -  HK)-1  H  A(r,{i})  (Ir  -  HK) 1  GU 

=  J  [Ir  +  (Ir  -  KH)-1  KH]  A(r,{i})  (Ir  -  HK)-1  GU 

=  J  (Ir  -  KH)-1  [Ir  -  KH  +  KH]  A(r,{i})  (Ir  -  HK)'1  GU 

=  J  (Ir  -  KH)-1  A(r, { i } )  (Ir  -  HK)-1  GU.  (6.3) 

Since  (9Y/9ki)k;  =  9Y/9lnlkil  and  K  A(r,{i})  =  kj  A(r,{i}),  it  follows  from  (6.3)  that 

9Y/3lnlkil  =  J  (Ir  -  KH)-1  K  A(r,(i})  (Ir  -  HK)-1  GU.  (6.4) 

Consider  taking  a  copy  of  the  system  in  Figure  6.1  with  zero  inputs,  inserting  a  summer  into  it, 
and  connecting  it  to  the  system  in  Figure  6.1  as  shown  in  Figure  6.2. 


Figure  6.2.  The  system  in  Figure  6.1  connected  to  a  duplicate  of  itself. 


Then  R<i  =  HKEd  +  V  =  HKRd  +  V;  hence,  Rd  =  (Ir  -  HK)'1  V  and 

Yd  =  JKRd  =  JK  (Ir  -  HK)-1  V  =  J  (Ir  -  KH)1  KV. 

From  Vfz)  =  A(r, { i) )  R(z)  and  (6.1), 

Yd  =  J  (Ir  -  KH)-1  K  A(r, { i } )  R  =  J  (Ir  -  KH)*1  K  A(r,{i))  (Ir  -  HK)-1  GU.  (6.5) 
Comparing  (6.5)  with  (6.4), 

Yd(z)  =9Y(z)/5lnikil.  (6.6) 

Since  the  operations  of  Z-transforming  and  differentiation  with  respect  to  a  parameter  can  be 
interchanged,  Yd(mT)  =  9Y(mT)/9lnlki!.  As  additional  notation,  let  Yd(z.S)  be  the  output  of  the 


duplicate  system  in  Figure  6.2  when  V(z)  =  A(r,S)  R(z).  In  this  notation,  Yd(z)  in  Figure  6.2  is 
Yd(z,{i}).  By  linearity, 

Yd(mT,S)=  X  3Y(mT)/3lnlkjl.  (6.7) 

ieS 

It  may  be  desired  to  set  kj  =  k  for  i  e  S  £  {l,2,...,r)  and  determine  3Y(mT)/3lnlkl  to 
examine  the  effect  of  changing  the  parameters  ki,  k  e  S,  in  unison.  Treating  each  kj  as  a  function 
of  k, 

f  k  ,  i  e  S 

ki(k)  = 

l  c;  ,  i  e  S 

where  C;  are  constants  independent  of  k,  and  using  the  chain  rule  for  differentiation, 

r 

3Y(mT)/3lnlkl  =  X  [@Y(mT)/3lnlkil)  OlnlkLl/ainlkl)] 

i=t 

=  X  aY(mT)/9lnlk,l. 
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From  (6.7), 

3Y(mT)/3lnlkl  =  Yd(mT,S).  (6.8) 

In  particular,  if  K  =  diag[k,k k]  then  S  =  { 1 ,2,...,r)  and  3Y(mT)/3lnlkl  is  found  by  setting 

V(mT)  =  R(mT).  By  a  similar  argument,  it  follows  that  if  the  parameter  k  appears  at  several 
locations  in  the  original  system,  BY (mT)/3lnlkl  will  be  the  output  of  the  duplicate  system  if  the 
signal  entering  the  gain  block  k  is  extracted  from  the  original  system  at  each  occurrence  of  k  and 
summed  with  the  signal  entering  the  gain  block  k  at  each  corresponding  occurrence  of  k  in  the 
duplicate  system.  Note  that  if  ki  and  k2  are  two  parameters  which  are  to  be  varied  independently, 
the  original  system  and  two  copies  of  it  are  required  to  determine  3Y(mT)/3lnlki  I  and 
3Y(mT)/3lnlk2l  simultaneously. 

The  quantity  3Y(mT)/3lnlkil  is  referred  to  as  a  semirelative  output  sensitivity 
function  in  the  literature  (Frank,  1978).  With  reference  to  Figure  6.2,  the  output  sensitivity 
function,  3Y(mT)/aik;l,  may  be  obtained  by  summing  V(mT)  with  the  output  of  the  gain  block  K 
in  the  duplicate  system. 

6.2  Interpretation  of  the  Time-invariant  Result 

Consider  a  multirate  discrete-time  system  with  input  u  and  output  y  which  contains  a 
discrete-time  subsystem  that  is  a  SISO  periodic  gain  g(nTj)  with  input  us,  output  ys,  and  period 
T:  =  qTj.  Thus,  this  gain  sweeps  cyclically  through  a  set  of  values  (gi . gq}  and  the 


subsystem  can  be  described  by 

ys(nT!)  =  g(nTi)  u,(nTj)  =  g((n  „,odq)+i)  uf(nTO. 

Let  T/T’  =  si,  where  T  is  the  period  of  the  entire  system  and  S2  e  N.  From  (2.6),  the  T-expanded 
representation  of  this  subsystem  is  Ys(mT)  =  K  Us(mT),  where  with  G  =  diag[gi,g2,...,gq], 
K  =  diag[G,G . G];  G  appears  S2  times  along  the  diagonal  of  K. 

In  light  of  the  time-invariant  nature  of  T-expanded  representations  and  the  fact  that  the 
structure  of  the  block  diagram  of  a  T-expanded  representation  is  identical  to  that  of  the  block 
diagram  of  the  original  lepresentation,  it  is  seen  that  the  T-expanded  representation  of  the  system 
under  consideration  here  is  amenable  to  the  application  of  the  variable  component  method  as 
described  in  the  preceding  section.  This  entails  passing  the  signals  entering  the  gain  K  in  the  block 
diagram  of  the  T-expanded  representation  of  the  system  through  a  matrix  gain  and  summing  the 
output  of  this  matrix  gain  with  the  input  of  the  gain  K  in  a  duplicate  of  the  T-expanded 
representation  of  the  system.  To  be  consistent  with  the  preceding  section,  let  r  =  qS2-  The  r 
parameters  composing  K  cannot  be  varied  independently  since  each  parameter  g;  appears  at  S2 
places  in  K.  Let  K  be  represented  as  K  =  diag[k},k2,...,kr],  and  for  each  i  6  { 1,2,.. .,q),  let  S;  be 
the  set  of  indices  of  the  kj's  which  correspond  to  the  positions  occupied  by  gp 

S,  =  (j  e  {1 . r}  :  kj  =  gi). 

From  (6.8),  it  follows  that  if  the  T-expanded  representation  of  the  original  system  (with  input 
L'(mT))  is  linked  by  the  matrix  gain  A(r,S,)  to  a  duplicate  of  itself  (with  zero  input)  as  described 
above,  then  the  output  of  the  duplicate  system  will  be  9Y(mT)/dlnlgil. 

A  system  in  real  time  whose  T-expanded  representation  is  that  proposed  above  to  produce 
5Y(mT)/3lnlgil  can  easily  be  found.  For  W  £  {l,2,...,q},  let  8(nTi,q,W)  be  the  Tiq-periodic 
discrete-time  function  of  n  given  by 


5(nTi,q,W) 


1  ,  if  (1  +  (n  mod  q))  €  W 
0  ,  otherwise 


Assume  T/qTj  =  S2  as  above;  from  (2.6),  a  SISO  discrete-time  periodic  gain  with  input  ug  and 
output  vg  described  by  yg(nTi)  =  8(nTi,q,{i})  ug(nTi)  is  seen  to  have  the  T-expanded 
representation  Yg(mT)  =  diag[A(q,{i)),...,A(q,{i))]  Ug(mT),  where  A(q, { i } )  appears  s^  times 
along  the  diagonal.  From  K  =  diag[G,...,G],  diag[A(q,{i)),...,A(q,{i})]  =  A(r,S;).  Thus,  if  the 


•  signal  entering  the  periodic  gain  g(nTi)  in  the  original  multirate  discrete-time  system  is  sampled  at 

!££  rate  1/Tj,  passed  through  the  periodic  gain  S(nT i ,q, { i } ),  and  summed  with  the  signal  entering  the 

V'  gain  g(nTi)  in  a  copy  of  the  original  system  (with  zero  inputs,)  the  T-expanded  representation  of 


n 


this  composite  system  coincides  with  that  described  previously  for  determining  3Y(mT)/8lnlg;l. 
Since  (3/3lnig,l)  can  be  distributed  across  the  components  of  the  vector  Y(mT),  5Y(mT)/01nlg;l  is 
the  T-expanded  version  of  0y/dln!g,l.  The  meaning  of  3y/9lnlgjl  stems  from  the  satisfaction  of 


ft 


-\  • 
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Assumption  2.1(f);  each  output  line  yj  of  the  system  is  sampled  and  held  at  some  rate  1/Tj  and  the 
corresponding  output  in  the  duplicate  system,  ydj(nTj),  gives  3yj(nTJ)/3lnlgiL 

The  following  procedure  summarizes  the  variable  component  method  for  multirate  discrete¬ 
time  systems.  This  procedure  is  more  general  than  the  results  discussed  in  this  section,  but  it 
follows  from  linearity  and  the  chain  rule  for  differentiation. 

Procedure  6.1: 

a.  Let  the  system  have  a  given  input  uo  and  resulting  output  yo.  Verify  that 
Assumption  2.1  is  met,  with  the  exception  that  periodic  single-rate  and  multirate 
discrete-time  subsystems  may  be  present. 

b.  Suppose  there  is  a  periodic  discrete-time  gain  g(nTj)  in  the  system  with 
period  qTi  and  having  values  (gi,...,gq),  as  described  above.  g(nTi)  may  be 
either  a  purely  discrete-time  gain  (such  as  a  parameter  in  a  computer  program,)  or 
the  discretized  version  of  a  continuous-time  gain  with  sample  and  hold  devices  at 
rate  1/Ti  on  its  input  and  output.  Let  W  £  { l,2,...,q),  and  for  all  i  e  W,  let  g;  be 
varied  in  unison  and  have  value  g.  The  parameters  gi,  i  g  W,  will  be  assumed 
independent  of  g. 

c.  Sample  the  signal  entering  g(nTi)  on  the  block  diagram  of  the  system 
with  rate  1/Ti.  Pass  this  sampled  signal  through  the  periodic  gain  5(nTi,q,W)  and 
sum  the  output  of  this  gain  with  the  signal  entering  g(nTj)  in  a  duplicate  of  the 
original  system. 

d.  Repeat  steps  (b)  and  (c),  using  the  same  duplicate  system,  at  each 
periodic  discrete-time  gain  in  the  system  that  has  parameters  with  value  g  if  these 
parameters  are  to  be  varied  in  unison  with  the  parameters  in  g(nTj).  Both  the  rate 
and  period  of  these  gains  may  be  different  from  location  to  location. 

e.  Simulate  or  implement  the  composite  system  resulting  from  the  above 
steps.  The  input  to  the  duplicate  system  should  be  set  to  zero,  and  the  initial 
conditions  on  both  the  original  and  duplicate  systems  must  be  zero.  The  output  of 
the  duplicate  system  is  then  9yo/01nlgl. 

f.  Add  to  the  composite  system  an  additional  duplicate  system  connected  to 
the  original  system  for  each  additional  parameter  that  is  to  be  varied  independently  if 
all  semirelative  output  sensitivity  functions  are  to  be  obtained  simultaneously. 

♦  ♦  ♦ 

The  variable  component  method  described  in  Procedure  6.1  possesses  properties  which 
simplify  the  mechanics  of  applying  Procedure  6.1.  With  reference  to  parts  (b)  and  (c)  of 
Procedure  6.1,  note  that  8(nTi,q,W)  =  1  only  at  the  times  when  g(nTj)  =  g.  This  is  seen  to  hold 
in  general;  at  each  discrete-time  gain  in  the  system,  values  are  passed  to  the  duplicate  system  only 


at  times  when  the  parameter  in  question  is  “active”  at  that  location.  In  particular,  if  g(nT[)  is  time 
invariant,  8(nTi,q,\V)  =  1.  Neither  the  T-expanded  representation  nor  the  value  of  T  is  needed  to 
apply  Procedure  6.1.  As  long  as  the  resulting  system  satisfies  Assumption  2.1  or  its  extensions, 
any  number  of  signals  in  the  system  may  be  considered  as  outputs,  and  these  outputs  may  be 
sampled  at  any  rate.  A  change  in  the  choice  of  outputs  or  their  rates  does  not  induce  any  change  in 
the  application  of  Procedure  6.1. 

As  with  the  variable  component  method  for  time-invariant  parameters,  Procedure  6.1  has 
the  property  that  an  additional  duplicate  system  is  required  for  each  independently  varied 
parameter.  If  each  value  of  a  periodic  gain  is  to  be  varied  independently,  the  number  of  duplicate 
systems  and  the  associated  computation  may  become  excessive.  Since  the  resulting  systems  will 
be  periodic,  it  is  in  general  insufficient  to  tune  the  response  of  the  system  to  a  fixed  input,  such  as  a 
step  at  time  t  =  0.  Instead,  it  may  be  necessary  to  tune  the  response  of  the  system  to  a  given  type 
of  input,  such  as  a  step,  applied  at  many  individual  times  during  the  period  of  the  system.  A 
system  with  input  uCnT^  exhibits  T/Ti  step  responses  which  are  in  general  distinct. 


6.3  Application  of  the  Variable  Component  Method 

The  semirelative  output  sensitivity  function  may  be  used  to  tune  the  time-domain 
performance  of  a  controlled  system  either  intuitively  or  through  the  solution  of  an  optimization 
problem.  The  variable  component  method  described  in  Procedure  6.1  provides  a  means  of 
obtaining  this  sensitivity  function  and  has  the  advantages  of  being  easy  to  apply  and  allowing  the 
tuning  of  periodic  gains  at  multiple  rates.  The  computational  burden  of  applying  Procedure  6.1 
may  be  substantial  if  the  number  of  parameters  to  be  varied  independently  is  large.  Conditions 
under  which  an  iterative  tuning  scheme  based  on  output  sensitivity  functions  involving 
minimization  of  a  cost  function  will  converge  or  result  in  a  stable  system  are  in  general  unknown  as 
of  the  present. 

Controllers  designed  with  the  aid  of  the  variable  component  method  offer  opportunities  for  i 

taking  full  advantage  of  multirate  sampling.  Since  Procedure  6.1  allows  outputs  to  be  sampled  at 
any  rate,  the  parameters  of  a  single-rate  controller  for  a  continuous-time  system  can  be  adjusted 
according  to  a  cost  based  on  samples  of  the  controlled  system’s  output  at  a  rate  which  is  faster  than 
the  rate  of  the  controller.  Any  accessible  signal  in  the  system  can  be  considered  as  an  output, 
enabling  the  parameter  changes  to  reflect  control  effort.  The  controllers  designed  in  Chapter  5 
share  the  property  that  they  are  not  directly  implementable  in  a  parallel  form  which  capitalizes  on 
the  multirate  nature  of  their  inputs  and  outputs.  An  application  of  the  variable  component  method  j 

might  involve  tuning  multiple  feedback  loops  containing  single-rate  controllers  of  a  fixed,  simple 
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structure  (such  as  PID,)  whose  rates  are  commensurate  with  the  time  scales  present  in  the  system. 
Such  controllers  are  readily  implemented  in  hardware  as  independent  parallel  processing  units  or  in 
software  as  concurrent  tasks. 

Example  6.1: 

In  Example  2.7,  assume  that  it  is  desired  to  tune  the  system  by  varying  gi 
and  g2  independently.  To  gain  a  more  accurate  picture  of  the  behavior  of  the  double 
integrator,  define  an  additional  output  of  the  system  as  yo,  the  output  of  the  double 
integrator  sampled  with  period  T4  =  0.1  sec.  As  y  is  a  sampled  version  of  yo,  it 
will  henceforth  be  ignored.  Application  of  Procedure  6.1  yields  the  composite 
system  shown  in  Figure  6.3,  where  yj  =  9yo/5lnlgil  =  (5yo/dgi)gi  and 
y  2  =  3yo/3Inlg2l. 


Figure  6.3.  Composite  system  resulting  from  Procedure  6.1. 


Denote  the  value  of  the  output  yo  at  time  t  =  mT4  with  parameter  values  gi  and  g2 
by  yo(mT4,g1,g2).  yi(mT4,g1,g2)  and  y2(mT4,glfg2)  are  defined  similarly.  At 
this  point,  a  number  of  options  are  available  for  tuning  the  system.  With  the 
intention  of  achieving  reasonable  settling  time,  overshoot,  and  rise  time,  consider 
choosing  gi  and  g2  to  minimize  the  cost 

200 

J  (gi.g2)  =  X  (mT4  [yo(mT4,gi,g2)  -  u0(mT4)])  2.  (6.9) 

m=0 

yo(mT4,gi+Agi,g2+Ag2)  can  be  approximated  for  small  Agi  and  Ag2  by 

y0(mT4,gi,g2)  +  [3yo(mT4,gi,g2)/9gi]  Agi  +  [3yo(mT4,g1,g2)/9g2]  Ag2, 
which  can  be  written  as 

yoCmTa.g!^  +  (Agi/gi)[yi(mT4,g1,g2)]  +  (Ag2/g2)[y2(mT4,g1,g2)].  (6.10) 

Substituting  (6.10)  into  (6.9)  gives  an  approximation  of  J(gi+Agi,g2+Ag2)  which 
is  quadratic  in  Agi  and  Ag2  and  can  be  minimized  analytically  in  terms  of  uo(mT4), 
yo(mT4,gi,g2),  yi(mT4,g!,g2),  and  y2(mT4,g1,g2)  which  result  from  simulation  of 
the  system  in  Figure  6.3  for  fixed  uo,  gi,  and  g2. 

Suppose  the  initial  parameter  values  are  gi(0)  =  -1.5  and  g2(0)  =  1  (these 
parameter  values  stabilize  the  system)  and  the  input  uo  is  a  3  sec  pulse  at  time  t  =  0 
of  unit  height.  The  response  of  the  composite  system  with  gi( 0)  and  g2(0)  to  uo  is 
shown  in  Figure  6.4.  Minimizing  J(gi(0)+Agi(0),g2(0)+Ag2(0))  gives  as  optimal 
values  Agi(0)/gi(0)  =  0.2633  and  Ag2(0)/g2(0)  =  0.4645  and  hence 

gi(l)  =  gi(0)  +  Agl(0)  =  -1.895  and  g2(l)  =  g2(0)  +  Ag2(0)  =  1.465. 

As  shown  in  Figure  6.5,  this  choice  of  parameter  change  produces  a  weighted 
combination  of  yi  and  y2  quite  suited  to  reducing  the  magnitude  of  the  oscillation  in 
yo-  Figure  6.6  shows  yo(mT4,gi(l),g2(l)),  which  is  seen  to  have  damping  much 
improved  over  yo(mT4,g1(0),g2(0)). 

After  repeating  the  process  of  minimizing  J(gi+Agi,g2+Ag2)  a  total  of 
seven  more  times,  the  values  of  gi  and  g2  for  all  practical  purposes  converge  to 
gi*  =  -2.714  and  g2*  =  2.008.  The  response  of  the  system  to  uo(nTi)  with  these 
parameter  values  is  shown  in  Figure  6.7.  The  state  transition  matrix  for  this  system 
is  unchanged  from  that  of  Example  2.7(i).  Substituting  gi*  and  g2*  into  this  state 
transition  matrix  gives  a  matrix  with  eigenvalues  -0.443  and  0.021.  The  final  tuned 
system  is  thus  asymptotically  stable.  Since  T/Ti  =  2,  this  system  exhibits  two  step 
responses.  Figure  6.8  shows  the  response  of  the  system  with  parameter  values  gi* 
and  g2*  to  uo(nTi-Ti).  If  the  response  in  Figure  6.8  is  deemed  undesirable,  one 


might  consider  repeating  the  tuning  process  by  minimizing  J(gi+Agi,g2+Ag2)  for 
both  the  input  uo(nTi)  and  the  input  uo(nTi-Ti)  at  each  stage,  requiring  two 
simulation  runs  for  each  iteration. 

♦  ♦  ♦ 


time  (sec) 


Figure  6.4.  Response  with  gj(0)  =  -1.5,  g2(0)  =  1.0. 


time  (sec) 


Figure  6.5.  Response  with  gi(0)  =  -1.5,  g2(0)  =  1.0: 

Ayo  =  (Agi(0)/gi(0))yi  +  (Ag2(0)/g2(0))y2. 
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time  (sec) 


Figure  6.6.  Response  with  gi(l)  =  -1.895,  g2(l)  =  1.465. 


time  (sec) 


Figure  6.7.  Response  with  gi*  =  -2.714,  g2*  =  2.008. 


time  (sec) 


Figure  6.8.  Response  to  uo(nTj-Ti)  with  gi*  =  -2.714,  g2*  =  2.008. 


CHAPTER  7 


CONCLUSIONS 


The  results  detailed  in  this  work  demonstrate  the  usefulness  of  the  T-expanded 
representation  of  a  multirate  system.  A  detailed,  systematic  procedure  was  presented  which 
enables  the  T-expanded  representation  of  a  multirate  system  to  be  calculated  using  any  linear 
systems  software  package  capable  of  matrix  operations  and  single-rate  discretization.  A  certain 
amount  of  information  is  clearly  lost  by  expressing  a  multirate  system  in  T-expanded  form,  most 
notably  the  interperiod  state  values.  However,  it  was  shown  that  the  t-varying  representation 
inherits  the  following  important  qualities  from  its  corresponding  T-expanded  representation: 
stability,  controllability,  rcconstructibility,  stabilizability,  and  detectability.  Although  this  argument 
was  carried  out  for  the  T-varying  representation  presented  here  (or  more  generally,  the  M-varying 
representation,)  the  process  of  defining  corresponding  representations  and  proving  the  theorems  in 
Chapter  3  and  Chapter  4  could  be  carried  out  for  almost  any  periodic  representation  of  the  class  of 
multirate  systems  treated  here.  The  multirate  controller  designs  in  Chapter  5  and  the  extension  of 
the  variable  component  method  presented  in  Chapter  6  demonstrate  that  the  T-expanded 
representation  and  the  notation  and  concepts  associated  with  it  find  practical  and  theoretical  use  in 
extending  time-invariant  results  to  multirate  systems.  Many  extensions  to  the  class  of  multirate 
systems  for  which  T-expanded  representations  were  defined  are  possible.  For  example, 
subsystems  employing  cyclic  sampling  (irregular  sampling  patterns  which  repeat  periodically  in 
time)  could  be  treated  by  altering  certain  connection  matrices  to  reflect  the  flow  of  information 
between  the  cyclic  samplers  and  the  other  samplers  in  the  system. 


APPENDIX 


SELECTED  PROOFS 


Proof  of  Theorem  2.1 

Lemma  A.l: 

If  the  sampling  periods  {Ti,...,Tn}  are  normalized  with  respect  to  TV, 
k€  {1,...,N},  then  GCD  {py  :  i  =  =  1. 

Proof:  (by  contradiction) 

Without  loss  of  generality,  suppose  that  normalization  is  performed  with 
respect  to  Tj.  Observe  that  rn  =  qn  =  1,  since  Ti  =  (l)Ti.  The  case  N  =  1  is 
trivial,  so  consider  only  N  >  1. 

GCD  {Ri/rn,...,Ri/riN}  =  1  because  Rj  =  LCM  {r^  :  i  =  1 . N}. 

Assume  that 

GCD  {pii  :  i  =  1,...,N)  =  GCD  {Ri,  Riqi2/ri2,-.„  RiqiN/riN}  =  n  *  1. 

For  this  to  be  true,  n  must  divide  Ri.  From  GCD  {Ri/rn,...,Ri/riN)  =  1,  there 
must  be  some  j  6  (2,3,...,N)  such  that  n  is  UQI  a  factor  of  Ri/rjj,  and  (because  n 
divides  Ri)  n  is  a  factor  of  rij.  For  GCD  {Riqii/rii :  i  =  1,...,N}  =  n  to  be  true,  n 
must  be  a  factor  of  Riqij/rij,  so  n  is  a  factor  of  qij.  This  is  a  contradiction, 
because  n  *  1  and  each  rii/q^  had  common  factors  removed. 

♦  ♦  ♦ 

Theorem  2.1: 

If  the  sampling  periods  are  normalized  with  respect  to  Tk  (giving  Tj  =  piuTk, 

i  =  1,...,N)  and  again  with  respect  to  Tj  (giving  T;  =  pjjXj,  i  =  1 . N,)  then  x k  =  tj 

and  pki  =  pji,  i  =  1,...,N. 

From  PkiTk  =  T;  =  pjiTj,  i  =  1,...,N,  it  follows  that 
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II  =  flu.  _  £k2  _  _  PkN 

Pjl  Pj2  PjN 

Thus,  ij/tk  is  a  rational  number.  Let  x/Xk  =  n/m,  where  n,m  e  N  and  n  and  m 
have  no  common  factor  other  than  1.  From 

_  rnpki  _  _  mpkN 

Pjl  PjN 

n  is  a  common  divisor  of  {pkjm  :  i  =  Since  n  and  m  have  no  common 

factors  other  than  1,  n  is  a  common  divisor  of  {pfc  :  i  =  By  Lemma  A.l, 

n  <  1.  Thus,  n  =  1.  Replacing  n  with  m  and  k  with  j  in  the  above  argument  yields 
m  =  1.  The  statement  of  the  theorem  then  follows. 

♦  ♦  ♦ 


Proof  of  Theorem  2.2 

Lemma  A.2: 

Given  a  set  of  sampling  periods  {Ti,...,Tn)  which  are  expressed  as 
T,  =  s;5,  where  Se  R,  s,  e  N  for  each  i,  and  LCM  { Si  :  i  =  1,...,N}  =  S,  the 
least  common  period  of  {Ti,...,Tn)  is  S8. 

Proof:  (by  contradiction) 

Assume  LCP  {Ti,...,Tn}  =  To  <  S8,  and  write  To  =  §8,  where  §  €  X  and 
§  <  S.  Then  each  sampler  will  sample  an  integral  number  of  times  during  the  time 
§8; 

§8  §8  ft 

T-  =  ”7  =  reNfori  =  1 . N- 

1 1  s,8 

This  implies  that$  >  LCM  {s, :  i  =  1,...,N)  =  S,  which  is  a  contradiction. 

♦  ♦  ♦ 


Theorem  2.2: 

T  =  x  (LCM  { pi  :  i  =  1,...,N})  =  Px  is  the  shortest  length  of  time  over 
which  a  multirate  system  with  sampling  periods  {Ti,...,Tn}  is  periodic. 

Proof:  A  trivial  application  of  Lemma  A.2. 

♦  ♦  ♦ 
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Proof  of  Theorem  2.3 


Theorem  2.3: 

Let  the  LCP  associated  with  {Ti,...,Tk}  be  T(k)  and  that  associated  with 
{T1,...,Tk,Tk+1}  be  T(k+1).  Then  T(k+1)/T(k)  e  N. 

Proof: 

Normalization  of  {Ti,...,Tk,Tk+i}  yields  T;  =  p^,  p4  €  N  for  i  =  l,...,k,k+l. 

By  Lemma  A.2,  T(k)  =  x(LCM  {pi,...,pk})  and  T(k+1)  =  x(LCM  {pi,...,pk+i }). 

Since  LCM  {pi,...,pk+i)  =  LCM  {LCM  {pi,...,pk},  Pk+i), 

T(k+1)  _  X(LCM  (pi,...,pk+i))  LCM  {LCM  (pi . pk ) .  Pk^i  I  N 

T(k>  x(LCM  {pi,...,pk})  LCM  {pi,...,Pk) 

♦  ♦  ♦ 
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